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ABSTRACT 

This  thesis  is  concerned  with  the  problem  of  computer  oriented 
objective  formulation  for  geographic  site  location  and  the  analysis  required 
to  locate  one  or  more  optimal  sites  based  upon  that  objective. 

For  most  practical  purposes,  a  geographic  region  is  a  two-dimen- 
sional Euclidian  space.   A  site  within  a  region  can  either  be  a  single 
point,  a  continuous  curve  or  path  between  two  points,  or  any  continuous 
subset  of  points  that  forms  an  area  within  the  region.   This  paper  will 
deal  primarily  with  the  location  of  areas.   It  will  be  shown  that  the 
problem  of  locating  an  area  can  subsume  the  problems  of  path  and  point 
location. 

Development  of  this  thesis  is  in  three  parts.   First,  another 
definition  of  areas,  paths,  and  points  is  proposed  in  terms  of  contiguous 
geographic  parcels.   This  definition  is  formalized  using  linear  graph  theory 
and  a  graph  model  of  the  parcel  configuration  under  study.   The  division  of 
geographic  regions  into  parcels  is  compatible  with  the  ways  that  gener- 
alized geographic  based  information  retrieval  systems  store  data.   Second, 
integer  programming  formulations  of  the  site  location  problem  are  developed 
which  are  known  to  be  solvable  and  will  provide  an  optimal  solution.   The 
results  of  several  tests  using  both  a  cutting  plane  algorithm  and  an  im- 
plicit enumeration  algorithm  are  discussed.   Third,  language  constructs 
are  implemented  within  an  existing  geographic  information  system  in  order 
to  generate  coefficients  for  objective  functions  and  constraints.   In  order 
to  efficiently  test  the  language  constructs  used  to  locate  optimal  sites, 
a  small  interactive  system  was  developed  employing  a  suboptimal  site 
location  heuristic. 
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CHAPTER  1 
INTRODUCTION 

There  must  be  an  infinite  number  of  motivations  for  man  to 
geographically  locate  sites.   Reasons  for  "finding  places"  are  encompassed 
within  the  entire  breath  of  human  social,  political,  environmental,  and 
economic  endeavors.   However,  man's  ability  to  find  a  "best  site",  what- 
ever the  purpose,  is  limited  by  the  quantity  and  quality  of  the  geographic 
based  information  available,  his  capabilities  at  manipulating  and  analyz- 
ing that  information,  and  his  grasp  and  formulation  of  an  objective  to  be 
optimized.   Computers  have  long  played  successful  roles  in  overcoming  the 
limitations  imposed  by  the  quantity  of  data.   Recently,  computer  analysis 
techniques  are  being  developed  to  deal  exclusively  with  spatial  data  - 
especially  for  environmental  modeling  problems  [l,2,3,^»5].   This  thesis 
deals  with  the  problem  of  computer  oriented  objective  formulation  for 
site  location  and  the  analysis  required  to  locate  optimal  sites  based 
upon  that  objective. 

For  most  practical  purposes,  a  geographic  region  is  a  two-dimen- 
sional Euclidian  space.  A  site  within  a  region  can  either  be  a  single  point, 
a  continuous  curve  or  path  between  two  points,  or  any  continuous  subset  of 
points  that  forms  an  area  within  the  region.   This  paper  will  deal  primarily 
with  the  location  of  areas.   It  will  be  shown  that  the  problem  of  locating 
an  area  can  subsume  the  problems  of  path  and  point  location. 

One  way  of  informally  stating  the  problem  of  site  location  is 
given  an  objective,  choose  areas,  paths,  and  points  within  a  geographic  re- 
gion which  satisfy  all  the  constraints  set  forth  for  the  sites  and  optimize 
the  objective.   This  description  is  clearly  related  to  the  language  used  in 
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the  formulation  of  mathematical  programming  problems.   Since  much  is  known 
about  the  solution  of  these  problems,  a  viable  approach  to  the  solution  of 
the  site  location  problem  would  be  to  transform  it  into  a  mathematical  pro- 
gramming problem. 

Development  of  this  thesis  is  in  three  parts.  First,  a  definition 
of  areas,  paths,  and  points  is  proposed  in  terms  of  contiguous  geographic 
parcels.   The  division  of  geographic  regions  into  parcels  is  compatible 
with  the  ways  that  generalized  geographic  based  information  systems  store 
data.   Second,  an  integer  programming  formulation  of  the  problem  is  developed 
which  is  known  to  be  solvable  and  will  provide  an  optimal  solution.  Third, 
the  problem  of  coefficient  formulation  for  creating  objective  functions  and 
constraints  will  be  explored. 


CHAPTER  2 
A  DEFINITION  OF  OPTIMAL  AREAS,  PATHS,  AND  POINTS 

2.1  Parcel  Configurations  for  Geographic  Regions 

Because  the  description  of  a  geographic  region  at  each  point  in 
the  continum  of  a  Euclidean  space  is  impossible,  several  approaches  have 
been  developed  to  summarize  geographies  [l,2,U].  The  most  popular  is  to 
segment  a  region  into  a  manageable  number  of  nonoverlapping  parcels  such 
that  their  union  covers  the  entire  region.  Data  is  then  collected  from  the 
geographic  area  represented  by  the  parcel  and  is  usually  stored  as  one  logi- 
cal record  within  a  geographic  information  system.   This  configuration  can 
be  further  categorized  into  two  main  types:   uniform  grid  and  generalized 
parcel  [k] .      Figure  1  depicts  these  structures.   Each  structure  presents 
advantages  and  problems  to  the  designers,  data  collectors,  and  users  of 
geographic  information  systems.   The  use  of  a  particular  structure  will  de- 
pend upon  the  application  for  which  a  system  is  intended  and  the  budgetary 
restrictions  imposed  upon  data  collection  and  system  development. 

Uniform  grid  systems  superimpose  a  uniform  rectangular  or  other 
regular  grid  over  a  region  [6,7].   Because  of  its  uniformity,  geographic 
information  systems  using  this  structure  are  the  easiest  to  construct.  Record 
addresses  can  be  calculated  directly  from  the  geographic  reference  which 
speeds  retrieval.   However,  data  acquisition  is  hindered  because  data  must 
be  collected  according  to  the  specific  configuration  of  a  uniform  grid.   Con- 
sequently, a  grid  size  for  these  systems  must  be  chosen  before  data  collection 
begins.   If  the  grid  is  very  small,  overhead  will  occur  in  keeping  track  of 
the  large  number  of  cells  needed  to  cover  a  region.   If  the  grid  is  large, 
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Figure  1.  Uniform  Grid  and  Generalized  Parcel  Configurations. 


5 
the  resolution  of  the  data  for  urban  regions,  where  landuse  information  is 

extensive,  might  be  too  low  to  be  useful. 

Generalized  parcel  systems  make  no  a  priori  assumptions  about  the 
shape  of  any  parcels  in  the  system.   The  parcels  are  usually  described  as 
arbitrary  polygons.   In  the  past,  these  systems  have  been  developed  to  store 
land  ownership  data.   Presently,  they  are  beginning  to  be  explored  for  the 
storage  of  extensive  environmental  information  [8,9s10,ll].  The  multiple 
resolution  problem  described  previously  for  a  uniform  grid  system  can  be 
solved  using  generalized  parcels.   This  advantage  must  be  balanced  against 
the  cost  of  digitizers  and/or  scanners  which  are  required  for  interpreting 
parcel  boundaries.   Furthermore,  the  storage  and  retrieval  algorithms  become 
significantly  more  complex  because  record  addresses  cannot  be  directly  cal- 
culated from  the  geographic  reference. 

2 . 2  Forming  Areas,  Paths,  and  Points  From  Sets  of  Parcels:   A  Linear 
Graph  Model 

For  the  remainder  of  this  thesis  a  geographic  region  will  be  ab- 
stracted as  a  two-dimensional  collection  of  nonoverlapping  parcels.   Each 
parcel  will  have  a  fixed  area  greater  than  zero.   However,  the  shape  or  con- 
figuration of  the  parcels  will  have  no  bearing  on  the  algorithms  developed 
for  site  location. 

A  site  within  a  geographic  region  of  parcels  must  be  redefined  in 
terms  of  parcels.  Figure  2  displays  area,  path,  and  point  sites  in  a  gener- 
alized parcel  configuration.  A  precise  definition  will  be  presented  using 
linear  graph  theory.   In  general  terms,  an  area  is  any  collection  of  contigu- 
ous parcels;  a  path  between  two  parcels  is  any  contiguous  set  of  parcels 
which  form  a  "line  or  string"  between  the  two  parcels;  and  a  point  is  a  single 
parcel.   It  is  apparent  that  a  path  is  an  area  with  special  linear 


POINT 


PATH 
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Figure  2.   Area,  Path,  and  Point  Sites  in  a  Parcel  Configuration, 


features  while  a  point  is  just  a  degenerate  area.   These  properties  will 
"be  discussed  shortly  in  more  detail. 

Figure  3  presents  a  generalized  parcel  system.  A  linear  graph 
that  models  this  system  is  formed  by  placing  a  node  or  vertex  in  each  parcel 
and  a  link  or  edge  between  two  nodes  if  and  only  if  they  are  adjacent  parcels, 
The  integer  programming  solution  presented  in  chapter  3  will  depend  on  the 
linear  graphs  being  oriented  or  directed  so  that  each  link  is  transformed 
into  two  directed  links,  one  in  each  direction.   Areas,  paths,  and  points 
can  easily  be  described  in  terms  of  subgraphs  of  these  linear  graphs. 

A  point  is  any  single  node  of  the  linear  graph. 

A  path  from  one  parcel  to  another  in  a  generalized  parcel  system 
is  a  set  of  contiguous  parcels  that  are  connected  by  a  directed  path  in  the 
directed  linear  graph  model  of  the  parcel  system.   Figure  h   displays  direc- 
ted paths  in  a  parcel  system.  A  directed  path  in  an  oriented  graph  is  a 
sub-graph  which  has  an  equal  number  of  arcs  entering  and  leaving  each  node 
with  the  exception  of  the  first  node  which  has  an  extra  arc  leaving  and  the 
last  node  which  has  an  extra  arc  entering.   This  definition  has  been  used 
to  establish  integer  programming  formulations  of  the  "shortest  path  problem" 
[12]. 

It  is  important  to  note  that  an  optimal  path  between  two  parcels 
will  never  be  modeled  by  a  directed  path  that  contains  a  directed  circuit 
or  loop.   Figure  h   shows  that  a  loop  contains  parcels  that  are  unnecessary 
for  the  completion  of  a  path.   In  this  Thesis,  the  term  directed  path  will 
only  apply  to  those  containing  no  directed  circuits. 

A  subset  of  parcels  form  a  single  area  within  the  study  region  if 
and  only  if  the  nondirected  subgraph  associated  with  these  parcels  is  con- 
nected. This  subgraph  is  formed  by  deleting  all  the  nodes  and  incident  arcs 


Figure  3.   The  Linear  Graph  Model  of  a  Parcel  Configuration, 
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10 
to  these  nodes  which  do  not  belong  to  the  subset  of  parcels  as  shown  in 
Figure  5«   From  linear  graph  theory,  it  is  known  that  there  exists  at  least 
one  tree  in  any  connected  linear  graph.   Also,  if  a  graph  has  n  nodes,  n-1 
links,  and  contains  no  circuits  then  it  is  a  tree  and,  therefore,  connected 

[13]. 

A  directed  tree  is  defined  to  be  a  directed  connected  linear  graph 
of  k  nodes,  k-1  links  with  no  more  than  one  link  directed  into  any  node,  and 
no  directed  circuits.   It  is  easy  to  show  there  must  exist  at  least  one 
directed  tree  in  the  directed  transformed  graph  that  corresponds  to  every 
tree  in  a  nondirected  graph.   Choose  any  node  in  the  nondirected  graph  to 
act  as  a  root  for  a  tree  and  find  the  corresponding  node  in  the  directed 
graph.   Find  all  the  nodes  which  are  daughters  of  the  root  in  the  nondirected 
graph  and  locate  all  the  corresponding  nodes  in  the  directed  graph.   By  the 
construction  of  the  directed  graph,  there  will  be  links  directed  from  the 
root  to  each  of  the  daughter  nodes.   By  choosing  each  daughter  in  turn  as  a 
new  root  and  finding  their  daughters  we  can  repeat  the  matching  process  until 
a  directed  tree  is  found  for  each  nondirected  tree.   Consequently,  if  a  space 
is  created  that  contains  all  directed  tree  subgraphs  of  the  directed  trans- 
formation of  a  graph  then  the  space  represents  all  possible  areas  or  subsets 
of  contiguous  parcels.  A  directed  tree  subgraph  of  an  area  is  shown  in 
Figure  6. 

It  is  important  to  note  that  both  a  directed  path  (with  no  directed 
circuits)  and  a  single  node  are  also  directed  trees.  Because  of  this  proper- 
ty, a  linear  graph  integer  programming  formulation  of  the  area  location  prob- 
lem also  subsumes  the  path  and  point  problems. 
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Figure  5.   The  Connected  Subgraph  of  a  Subset  of  Contiguous  Parcels 
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2 . 3  Parcel  Weighting,  Objective  Formulation,  and  Opt imal  Site  Definition 

A  mechanism  must  be  provided  which  establishes  the  desirability  of 
one  site  over  another.   One  way  to  accomplish  this  is  to  numerically  weight 
each  parcel  in  the  region  by  a  positive  number  with  respect  to  the  contribu- 
tion, favorable  or  adverse,  which  a  parcel  would  add  should  it  be  chosen 
as  part  of  the  site.   The  value  of  a  site  will  simply  be  the  sum  of  all  the 
weights  of  the  parcels  that  make  up  the  site.   Methods  for  using  geographic 
based  information  systems  to  establish  these  weights  will  be  developed  in 
Chapter  k. 

Let  us  assume  that  meaningful  weights  can  be  established  from  the 
data  stored  for  the  parcels.   Then  the  objective  is  to  determine  the  site  of 
parcels  that  minimize  or  maximize  the  value  of  the  site.   If  the  weights 
represent  suitabilities  then  a  maximum  valued  site  will  be  sought.   However, 
if  the  weights  describe  costs  then  a  minimum  valued  site  will  be  best.   Fig- 
ure 7  shows  a  weighted  parcel  system  and  some  optimally  allocated  sites. 
Therefore,  given  a  parcel  configuration  for  a  region  and  a  weight  for  each 
parcel,  the  definition  for  an  optimal  site  becomes  "Maximize  or  minimize  the 
sum  of  the  weights  of  the  parcels  that  form  the  site  subject  to  the  con- 
straints that  l)  the  site  is  an  area,  path,  or  point  of  parcels  and  2)  that 
all  other  requirements  for  the  site  are  satisfied".   "Other  requirements" 
are  any  conditions  that  must  be  met  before  a  feasible  collection  of  parcels 
can  be  an  optimal  solution.   Conditions  might  be  prescribed  on  the  size  of 
an  area,  the  beginning  and  end  points  of  a  path,  or  any  parameter  of  the 
data  associated  with  the  site. 
2.U   Simultaneous  Multiple  Site  Location 

It  is  possible  to  encounter  situations  where  two  or  more  sites 
must  be  located  simultaneously  such  that  locational  interdependences  must 
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Figure  7-   Optimal  Area,  Path,  and  Point  Sites  in  a  Parcel  Configuration. 


Ik 

be  taken  into  account.  Consider  the  problem  of  locating  an  industrial  park 
and  a  large  residential  subdivision  to  house  the  park's  new  employees.  Lo- 
cating each  site  independently  would  not  take  into  consideration  the  undesir- 
ability  of  having  housing  too  close  to  noise  and  air  pollution  sources.  On 
the  other  hand,  excessive  transportation  time  and  cost  would  be  incurred  if 
housing  and  employment  were  distantly  separated. 

The  above  situation  as  well  as  similar  ones  can  be  handled  by  con- 
sidering extensions  to  the  site  location  objective.   Not  only  must  the 
weights  of  each  parcel  for  each  site  be  expressed  in  the  objective  but  also 
a  weighting  function  of  distance  between  parcels  in  dii'Serent  sites  for  each 
pair  of  sites  must  be  included. 

The  next  chapter  will  establish  the  integer  programming  activity 
variables  and  constraints  required  to  develop  solvable  formulations  of  the 
single  and  multiple  site  location  problems. 
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CHAPTER  3 
INTEGER  PROGRAMMING  SOLUTIONS 

The  site  location  algorithm  exploits  the  directed  linear  graph 
model  of  a  geographic  parcel  configuration.  A  set  of  linear  constraints  is 
developed  with  reference  to  the  graph  which  creates  a  space  of  feasible  in- 
teger points.   Each  point  describes  a  directed  tree  subgraph  that  represents 
a  possible  site.   These  constraints  are  then  inserted  into  an  integer  pro- 
gramming framework  where  the  objective  is  to  optimize  a  linear  function 
whose  coefficients  are  the  weights  of  the  corresponding  parcels.   The  solu- 
tion to  this  integer  programming  formulation  will  be  the  parcels  of  the 
sites,  represented  by  integer  points  in  the  feasible  space,  that  minimizes 
or  maximizes  the  objective  function. 
3 . 1  An_  All-Integer  Formulation 
3.1.1  Formulation  of  the  Initial  Tableau 

It  was  shown  in  section  2.2  that  a  directed  tree  is  defined  by  the 
following  properties: 

(1)  the  directed  graph  has  k  nodes  and  k-1  links; 

(2)  no  more  than  one  link  is  directed  toward  any  node; 

(3)  no  directed  circuits  are  formed. 

The  linear  space  of  all  possible  directed  trees  will  be  established  by  form- 
ing constraint  equations  that  mirror  the  above  properties. 

In  order  to  establish  the  constraint  equations,  integer  programming 
"activity"  variables  are  attached  to  the  directed  graph  of  a  parcel  system. 
Figure  8  shows  a  portion  of  a  directed  linear  graph  and  its  associated 
activity  variables.   For  each  node  i   (representing  parcel  i),  i=l,2,...,n, 
where  n  is  the  number  of  parcels  or  nodes,  let  N.  be  the  activity  of  including 
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a  .  ,  w . ,  K  .  ,  U  . 
J    J    J    J 


Figure  8.   The  Integer  Programming  Activity  Variables. 


IT 
node  i  in  a  feasible  tree  subgraph.   Let  L. .  be  the  activity  of  including 
the  link  directed  from  node  i  to  node  j   in  a  feasible  tree  subgraph. 
Also,  let  U.  be  a  utility  variable  for  node  i. 

Let, 

_  I  1,  if  node  i  is  in  the  solution  tree 
i    \0,  otherwise 

Si,  if  the  link  from  i  to  j  is  in  the  solution  tree 
ij   |0,  otherwise 

and 

U.=  a  nonnegative  integer 

All  integer  programming  variables  are  implicitly  nonnegative. 
Each  node   i  also  has  a  fixed  positive  weight  w.  and  area  a.. 

The  requirement  that  for  k  nodes  there  can  only  be  k-1  links  or, 
equivalently,  that  there  must  be  one  less  link  than  the  number  of  nodes  is 
satisfied  by 


I     N.-     Z  L. .  =  1  (1) 

i=l   X  all  i,j  s.t.   1J 
L.  .  in  graph 

The  second  requirement  that  no  more  than  one  arc  can  be  oriented  toward  any 
node  is  satisfied  by 

Z       L.k  <  1,  k=l,2,...,n  (2) 

all  i  s.t. 

Lik  in  6raph 

Thirdly,  the  requirement  that  there  exist  no  directed  circuits  is  fulfilled 
hy 

U.-  U.+  nL. .  <  n-1,  for  all  L..  in  the  graph.  (3) 
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This  set  of  equations  can  be  shown  to  eliminate  all  subgraphs  with  directed 
circuits  while  not  eliminating  any  feasible  directed  tree  subgraphs.   It 
has  been  used  in  the  past  as  part  of  an  integer  programming  formulation  of 
the  traveling  salesmen  problem  [12].   To  show  that  all  directed  circuits 
are  eliminated,  consider  any  directed  circuit  possible  through  k  nodes. 
This  requires  that  L.  =  1  for  k  of  the  links.  The  k  associated  con- 
straints  are: 

(0)    U.  -  U.  +  nL.   .   <  n-1; 

X0        Xl    V1!  " 


(1)    U.  -  U.  +  nL.   .   <  n-1; 
ix   i2    il9i2  " 


(k-2)  U.    -  U.    +  nL.     .    <  n-1; 
i      i       it    — 
k-2    k-1     k-2'  k-1 

(k-1)  U.   -  U.  +  nL.     .  <_  n-1. 
\-l        10    V-l^O 

It  is  well  known  that  for  any  solution  to  be  feasible,  it  must  also  satisfy 
the  constraint  formed  by  the  addition  of  any  other  constraints  in  the  model, 
Therefore,  the  circuit  formed  by  the  k  links  must  also  satisfy 


k-1 

I      (U.  -  U.  +  nL.   .  <_  n-1) 

j=0    j    j+1  mod  k      j'  j+1  mod  k 


k-1  k-1 

=>      Z      (U.    -  U.  +  nL.       .  )    <_       I      (n-1) 

j=0       ^        X j+1  mod  k  Wlmodk  j=0 


k-1  k-1  k-1 

=>  I      (U.  -  U.         )  +   Z  nL.   .  <_   Z   (n-1) 

j=0   Xj   1j+l  mod  k    j=0   Xj  iXj+l   mod  k  "   j=0 


Since  both  U.   and  -U.   for  j=0,l,2, . . . ,k-l  appears  in  the  leftmost  summa- 

1j       1J 
tion  above,  we  have 
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k-1 

I   (U.  -  U.         )  =  0. 
j=l   1j    j+1  mod  k 

Also,  since  L.  =1  for  j=0,l, . . . ,k-l,  the  above  inequality  becomes 

j '  j+1  mod  k 

k-1     k-1 

E  n  <   E   (n-1) 
j=0    "  j=0 

=>  kn  <_  k  ( n-1 ) 

However,  this  is  a  contradiction.   Thus,  any  directed  circuit  must  be  excluded, 

To  show  that  no  solution  tree  is  eliminated  by  constraint  set  (3), 
it  must  be  shown  that  there  exist  values  for  the  U. 's  that  satisfy  (3)  for 
any  given  solution  tree.   Consider  a  solution  tree  which  includes  k  nodes 
and  k-1  links.   The  root  of  the  tree  is  defined  to  be  position  1   and  each 
node  is  numbered  consecutively  from  left  to  right,  one  tree  level  at  a  time, 
and  with  the  lowest  levels  first.  Figure  9  displays  this  numbering  scheme. 
Let 

t  =  p,  the  p   position  in  the  tree  for  p=l,2,...,k 

ir 

and 

t  =  p,  arbitrarily  represent  the  n-k  nodes  not  in  the  solution  tree 
where  p=k+l,k+2, . . . ,n. 

Clearly,  there  exists  a  1-1  correspondence  between  the  t  's  and  the  N.'s. 
Also,  there  exists  a  1-1  correspondence  between  the  N. 's  and  U. 's  by  def- 
inition of  the  activity  variables.   Let 

U.=  p  where  U.<->  N.<->  t  =p,  i ,p=l,2, . . . ,n. 
i  *  i    i    p 

The  symbol  "<->"  denotes  correspondence. 

To  show  that  the  above  is  a  feasible  set  of  values  for  the  U. 's, 

l 

suppose  L   is  a  link  in  the  given  solution  tree.   Thus, 
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LEVEL 


LEVEL  2 


LEVEL  3 


LEVEL  4 


Figure  9-   A  Node  Numbering  Scheme  for  Solution  Trees. 
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and 


L  =  1 
rs 


N  =  N  =  1. 
r   s 

From  above  we  know  for  some  fixed  integers  a  and  b  between  1  and  k  that 

K  <_>  u  <->  t  and  N  <->  U  <->  t. 

r     r     a      s     s     b 

which  implies  by  definition  that 

U  =  a  and  U  =  b. 

r         s 

By  construction  of  the  tree  positions  t  and  t,  we  have 

a      b 

b>a  or  U  >  U  . 
s   r 

Thus,  the  constraint  which  includes  L   , 
'  rs 

U  -  U  +  nL   <  n-1, 
r   s    rs  — 

is  satisfied  since  U  >  U  and  L  =  1.   Now  suppose  L   is  a  link  that  is  not 

s   r      rs  rs 

in  the  given  solution  tree.   Thus, 

L  =  0 
rs 

causing  the  constraint  from  (3)  which  includes  L   to  become 
to  rs 

U  -  U  <  n-1 
r   s  — 

which  is  satisfied  because  U  <  n  and  U  >  1. 

r  —       s  — 

The  previous  constraint  sets  (l),  (2),  and  (3)  are  not  sufficient 

to  define  the  model.   Two  more  sets  of  equations  must  be  added.   First,  the 

K.  and  L. .  variables  must  be  constrained  to  0  or  1.   Equation  (2)  implies 

that  all  the  L..'s  are  less  than  1,  however,  the  following  must  be  added 

to  restrict  the  W.'s. 

l 

Ni  <  1,  i=l,2,...,n  (k) 
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Second,  the  L..  and  N.  variables  must  "be  forced  to  interact.   The  following 
equations  can  be  shown,  by  setting  up  the  exhaustive  table  of  possible  out- 
comes, to  force  both  N.=  1  and  N.=  1  when  L. .=  1  or  L..=  1.   In  addition, 

i         J  ij        Ji 

L. .=  L..=  1  is  a  directed  circuit  which  can  be  eliminated  as  necessary  to 
formulate  equation  set  (5)-   Table  1  displays  the  allowed  proper  combinations 
of  values  for  N. ,  N.,  L..,  and  L..  which  is  satisfied  by 

-2L.  .  -2L..  +  N.  +  N.  >  0,  for  all  L..,  L..  pairs  in  the  graph.     (5) 
ij    ji    l    J  ij   ji 

Equations  (l)  through  (5)  insure  that  the  optimal  solution  will 
be  a  directed  tree.   However,  additional  constraints  are  required  so  that 
the  directed  tree  will  represent  a  proper  area,  path,  or  point.   For  instance, 
the  size  of  an  area  must  be  defined  before  an  optimal  area  can  be  located. 
Recall  that  each  node   i  has  an  associated  area  a  •  •  The  constraints  on  the 
total  area  required  of  the  contiguous  parcels  are, 

n 

A„  <   Z  a.N.  <  A   ,  (6a) 

I  —  .    .      l  l  —  u 
1=1 

where  A  is  an  upper  bound,  and  A„  is  a  lower  bound, 
u  I 

A  path  site  location  tableau  is  created  by  adding  constraints  to 
sets  (l)  through  (5)  that  specify  the  end-points  of  the  path  and  that  force 
the  solution  tree  to  be  a  path.   The  end-points  can  be  specified  by  setting 
their  respective  node  variables  equal  to  one.   By  choosing  one  end-point  to 
be  the  root  of  the  tree  (path)  and  setting  the  sum  of  the  links  which  are 
oriented  away  from  the  root  equal  to  one,  we  will  necessarily  force  a  path 
between  the  end-points  when  minimizing.   It  should  be  observed  that  this 
summation  constraint  will  cause  the  root  to  be  forced  into  the  solution  since 
one  of  the  links  directed  from  it  will  be  in  the  solution.   This  is  guaran- 
teed by  equation  set  (5).   Thus,  setting  the  node  variable  for  the  root  equal 
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L.  . 

L.. 

N. 

l 

K. 
J 

-2L. .-2L..+  N.+  N 

Allowed 

0 

0 

0 

0 

0 

Yes 

0 

0 

0 

1 

1 

Yes 

0 

0 

1 

0 

1 

Yes 

0 

0 

1 

1 

2 

Yes 

0 

1 

0 

0 

-2 

No 

0 

1 

0 

1 

-1 

No 

0 

1 

1 

0 

-1 

No 

0 

1 

1 

1 

0 

Yes 

1 

0 

0 

0 

-2 

No 

1 

0 

0 

1 

-1 

No 

1 

0 

1 

0 

-1 

No 

1 

0 

1 

1 

0 

Yes 

1 

1 

0 

0 

-k 

No 

1 

1 

0 

1 

-3 

No 

1 

1 

1 

0 

-3 

No 

1 

1 

1 

1 

_______ 

-2 

No 

Table  1.   Allowed  Outcomes  of  N.,  N  ,  L   ,  and  L 

1      J       J- J  J  -*~ 


2k 
to  one  is  not  necessary.   Therefore,  only  the  following  constraints  must  be 

added  to  equations  (l)  through  (5)  to  locate  a  path  between  nodes  N.  and  N  . 

i      J 

Z  L  =  1 

all  k  s.t. 

L.,  in  graph 
lk 

and  ( 6b ) 

N.=  1 
J 

Notice  that  constraints  for  path  location  only  work  when  a  minimum  weighted 
path  is  sought.  Maximizing  the  objective  with  constraints  (l)  through  (5) 
and  (6b)  would  encompass  every  parcel  because  each  parcels'  weight  in  the 
region  would  add  to  the  value  of  the  site  and  there  may  be  no  restriction 
on  the  total  area  or  the  maximum  number  of  parcels  in  the  path.   In  order 
to  convert  from  maximization  to  minimization,  the  weights  w.  must  be  com- 
plemented by  making  the  largest  weights  smallest  and  the  smallest  weights 
largest.   This  can  be  done  by  finding  W  such  that 
W  =  max  {w. ,  i=l,2, . . . ,n} . 

The  complementary  weight,  w. ' ,  can  be  constructed  for  each  parcel  by 

w!  =  W  -  w.  +1,  i=l,2,...,n. 
l        l 

In  order  to  locate  a  single  point,  the  simple  constraint 

n 

Z     N.=  1  (6c) 

i=l  1 

could  be  added  to  equations  (l)  through  (5).   However,  equation  (6c)  alone 
would  be  sufficient  to  specify  the  integer  programming  formulation  since 
adjacency  relationships  play  no  role  in  point  location. 

To  conclude  the  model,  the  objective  function  necessary  to  find 
the  minimum  or  maximum  weighted  site  is 
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n 


Z  w.N.  (7) 

..11 

i=l 

where  the  optimal  site  is  formed  by  the  set  of  parcels 

{il  N.=  1}  . 
1   l 

3.1.2  A  Canonical  Formulation 

A  canonical  formulation  of  an  integer  programming  problem  is 

■+   ->  T 
minimize  c  x 

subject  to  the  constraints: 
Ax  £  b 

x  >^ 

x  e  Zn 
where  A  is  an  m  x  n  matrix,  c  e  R  ,beR  ,  R  is  the  set  of  all  reals,  and 
Z  is  the  set  of  all  integers  [lU]. 

If  an  integer  programming  problem  is  not  given  in  the  above  form, 
it  can  be  transformed  using  the  following  rules: 

(1)  Replace  all  equations  by  the  two  inequalities  <_  and  ^_; 

(2)  Multiply  inequalities  of  the  form  _>  by  -1; 

(3)  Multiply  the  objective  function  by  -1  when  maximizing. 
Table  2  displays  the  canonical  formulation  for  the  site  location  problem. 

An  example  will  be  shown  in  the  next  section  which  shows  how  the  constraints 
are  formulated  in  practice.   Formulas  for  the  number  of  variables  and  con- 
straints required  to  formulate  an  area  location  problem  in  a  p  x  p  rectangu- 
lar parcel  configuration  are: 

2 
number  of  variables  =  6p  -  Up 


and 


2 
number  of  constraints  =  8p  —  6p  +  1+ . 


OBJECTIVE  FUNCTION: 
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n 
MINIMIZE   Z  v.N. 
i-1   X  X 


SUBJECT  TO  THE  CONSTRAINTS: 


FOR  DIRECTED  TREES 


(1)   Z  N.  -      Z      L.  .  <_  1; 

i=l   1  all  i,j  s.t.   1J 

L. .  in  graph 
ij 

n 

Z       L.  .  -   E  N.  <  -1 

all  i,j  s.t.   1J   i=l  1  " 

L.  .  in  graph 


(2)     z        Lik  <_  1,  k=l,2,..  .,n 
all  i  s.t. 
L.,  in  graph 


(3)  U.-U.+nL..  <  n-1,  for  all  L. .  in  the  graph 


(k)     N.  <  1,  i=l,2,...,n 


(5)  2L..+2L..-N.  -N.  <  0,  for  all  L..,  L..  pairs  in  the  graph 


Table  2.  A  Canonical  FormixLation  of  the  Site  Location  Problem. 
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FOR  AN 

AREA 

n 

(6a) 

I     a.N. 
1=1   X  X 

n 

< 

A  ; 
u 

Z     (-a. 
1=1    X 

)N.  < 
l  — 

, 

FOR  A 

PATH  FROM  N.  TO  H. 

(6b) 

E        L.k<l; 
all  k  s.t. 

Lik  in  §raPh 

ii*    ♦      (-lik)--1; 

all  k  s.t. 
L   in  graph 

N.  <  1; 
J  - 

-N.  <  -1 
J  - 

FOR  A 

POINT 

n 

(6c) 

I      N.  <  1; 
i=l   1_ 

n 

I      ( -N .  )  < 
i=l    X 

1 

Table  2    (continued) 
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The  spare eness  of  the  initial  tableau  is  defined  as  the  number  of  nonzero 
coefficients  divided  by  the  total  number  of  coefficients  and  is  estimated  by 

37P2  -  32p _  -,  -2s 

sparceness  = *■ - — *- =  0(p   ). 

(6p^  -  Hp)(8p  -  6p  +  h) 
3.1.3  A  Mote  on  Additional  Constraints 

In  the  definition  of  an  optimal  site  given  in  section  2.3,  it  was 
stated  that  all  requirements  set  forth  site  must  be  met  in  addition  to  the 
area,  path,  or  point  constraints.  Requirements  can  be  handled  within  the 
integer  programming  framework  without  any  additional  analysis  if  they  can  be 
formulated  as  constraints  to  be  added  to  the  initial  tableau  of  constraints 
(l)  through  (6).   For  example,  a  particular  application  might  require  the 
optimal  site  to  be  located  where  the  total  population  of  the  parcels  must  be 
less  than  or  equal  to  P.   Let  p.  be  the  population  of  parcel   i.   The  con- 
straint is  written  by 

n 

I  p.N.  <  P. 
i=l  X  X  " 

3.1.U  An  Example 

Consider  the  weighted  parcel  configuration  shown  in  Figure  10. 
The  coordinate  triple  that  describes  each  parcel  is  comprised  of  the  parcel 
number,  area,  and  weight.   Let  the  example  problem  be  the  location  of  the 
minimum  weighted  area  such  that  the  total  size  of  the  site  is  between  h   and 
6  units. 

Part  of  the  integer  programming  initial  tableau  will  be  set  up  in 
terms  of  equations  (l)  through  (7)  of  section  3.1.2.   Other  rows  of  the 
tableau  are  required  by  Gomory's  cutting  plane  algorithm  which  was  used  in 
this  thesis  for  solving  all-integer  programming  problems  [12, lU],   The  term 


PARCEL   NO. 

AREA 

WEIGHT 

1 

2 

4 

2 

2 

2 

3 

1 

1 

4 

2 

2 

5 

3 

5 

6 

2 

5 

7 

1 

2 

8 

1 

3 

9 

1 

2 

10 
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"all-integer  programming"  refers  to  problems  in  which  all  the  coefficients 
in  the  initial  tableau  as  veil  as  the  feasible  values  of  the  activity  vari- 
ables are  integers. 

The  ALGOL  code  of  Gomory's  algorithm  expects  the  initial  tableau 
to  be  of  the  following  form  [15]: 

(1)  An  integer  matrix  A  =  [a.   ]  of  m+1  rows  numbered  0  to  m 
and  n  columns  numbered  1  to  n. 

(2)  The  elements  of  row  i,  i=l,2,...,  m-n+1 ,  represent  the 

constraint  a.  .x,  +  a.  ^x_  +  . . .  +  a .   n<a.   .  The  elements 
i,l  1    i,2  2         i,n-l  —  i,n 

of  the  last  column  of  A,  except  the  topmost,  are  seen  to  be 
the  right  hand  sides  of  the  constraint  inequalities. 

(3)  The  activity  variables  x  ,x  , . .  . ,  x    are  implied  by  column 
position. 

(h)     The  elements  of  row  0  represent  the  objective  function 

a.  ,  x  +  a_  „x  +  .  . .  +  a_   ..x  ..  which  is  to  be  minimized. 
0,1  1   0,2  2         0,n-l  n-1 

(5)  The  element  a0   must  be  initially  set  to  0  but  will  become 

0,n 

the  negative  of  the  value  of  the  objective  function  when  the 
algorithm  terminates. 

(6)  The  elements  of  row  i,  i=m-n+2, . . . ,n,  are  assigned  as  follows: 

f  -1,  if  j  =  i-m+n-1 
ij    1  0,  otherwise 

These  rows  represent  the  constraints 

-x.  <  0  or  x.  >  0,   i=l,2,...,  n-1, 
l  —        l  — 

which  are  not  implied  by  Gomory's  algorithm. 

(T)  The  algorithm  will  terminate  if  all  entries  in  the  last 

column,  except  the  topmost,  are  non-negative.   The  optimal 

solution  is  given  by  x.=  a.  ,    , .,   ,  i=l,2, . . .  ,n-l. 

i   i+m-n+l,n 
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(8)  All  the  tableau  columns,  except  the  last,  must  be  lexicogra- 
phically positive.   This  means  that  the  first  nonzero  element 
in  these  columns  must  be  positive. 
The  program  which  generates  initial  tableaus  and  includes  the 
Gomory  integer  programming  algorithm  is  given  in  Appendix  A.   In  order  to 
generate  the  initial  tableau,  a  tabular  description  of  the  directed  linear 
graph  model  of  the  parcel  configuration  must  be  provided.   Each  link  is 
assigned  an  identifying  number  in  the  table.   Figure  11  displays  the  linear 
graph  for  the  example  problem.   Table  3  shows  the  tabular  description  of  the 
linear  graph  model.   Each  line  of  the  table  is  converted  to  a  card  image  for 
input  to  the  program.   The  Table  is  converted  by  the  program  into  a  two- 
dimensional  array  ADJ  where 

(identification  number  of  the  link  from  node  I  to  J,  if  I 
ADJ[l,J]  =<is  adjacent  to  J 
jO,  otherwise 

Constraints  (rows  of  matrix  A)  are  generated  by  procedures  that 
refer  to  array  ADJ.   Each  procedure  corresponds  to  one  of  the  equation  sets 
(1)  through  (7). 

Figure  12  displays  a  listing  of  the  initial  tableau  and  final 
solution  for  the  example  problem.   Figure  13  shows  how  the  tableau  is  divided 
into  equation  types.   Note  that  the  node  variables  are  represented  by  the 
first  10  columns.   The  link  variables  are  located  in  columns  numbered  by 
adding  10  (the  number  of  nodes)  to  the  link  number  given  in  array  ADJ.  Fi- 
nally, the  utility  variables  correspond  (by  column  number)  to  the  corres- 
ponding node  variables  by  adding  36  (the  number  of  links  and  nodes)  to  the 
node  number. 

Row  38  represents  constraint 
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Figure  11.   The  Linear  Graph  Model  of  the  Example  Problem. 
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Figure  12.   The  Initial  Tableau  and  Final  Solution  of  the  Example  Problem. 
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Figure  13.   The  Equation  Set  Structure  of  the  Example  Initial  Tableau. 


37 

n         /    % 

E  U.  < ,  n  =  the  number  of  nodes, 

.  ,   1  —     2 
1=1 

which  is  incorporated  to  meet  the  requirement  that  the  initial  tableau  be 

lexicographically  positive.  Recall  that  the  proof  of  equation  set  (3)  in 

section  3.1.1  showed  that  it  was  sufficient  to  let  each  utility  variable  be 

assigned  a  different  integer  in  order  for  the  constraints  to  fulfill  the 

function  of  restricting  all  directed  circuits.  The  above  constraint  allows 

each  of  the  n  utility  variables  be  assigned  a  different  integer  from  1  to  n. 

This  results  in  the  summation  of  l+2+...+n  = as  an  upper  bound  on 

the  summation  of  the  utility  variables. 

3.1.5  Computational  Experience 

Computational  tests  were  run  on  small  uniform  grid  configurations 
in  order  to  demonstrate  the  correctness  of  the  integer  programming  formula- 
tion as  well  as  to  determine  the  computational  efficiency  of  the  approach. 
A  published  ALGOL  code  for  Gomory's  cutting  plane  algorithm  was  used  to  test 
the  formulation  [15] •   Testing  was  performed  on  a  Burroughs  B67OO  [l6].  Al- 
though the  code  is  a  very  inefficient  implementation  of  Gomory's  method  tak- 
ing no  advantage  of  the  special  structures  found  in  the  initial  tableau,  it 
does  provide  insite  into  computational  characteristics  of  the  model. 

Weights  were  arbitrarily  chosen  for  several  p  x  p  uniform  rectangu- 
lar grid  parcel  configurations.   Figure  lU  shows  the  average  results  of  sev- 
eral problems  of  varying  size.   The  computational  time  is  seen  to  be  directly 
proportion  to  both  the  size  of  entire  region  as  well  as  the  size  of  the  site 
within  the  region.   It  was  also  noted  that  the  times  were  larger  for  problems 
in  which  several  alternative  optimal  solutions  existed  for  a  site. 

Other  users  of  integer  programming  have  reported  that  the  ordering 
of  constraints  in  the  initial  tableau  has  an  effect  on  processing  time.   This 
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Figure  1*4 .   The  Average  Processing  Times  of  Several  p  x  p  Problems 
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was  explored  with  very  favorable  results  —  as  much  as  an  order  of  magnitude 

reduction  occurred.   The  constraint  order  shown  in  the  example  tableau  of 

figure  13  was  found  to  be  best  for  all  test  problems. 

Extensions  to  the  code  for  Gomory's  method  have  been  made  to 
incorporate  the  potentially  stronger  cuts  derived  by  Wilson  [IT]-   However, 
tests  using  the  extended  code  showed  no  appreciable  gains  or  losses  for  any 
of  the  test  problems . 

It  has  long  been  known  that  Gomory's  algorithm  for  solving  integer 
programming  formulations  has  demonstrated  excessive  computational  times. 
Unfortunately,  the  model  developed  in  this  thesis  is  no  exception.   It  is 
easy  to  see  that  large  scale  problems  of  hundreds  of  parcels  would  generate 
problems  of  thousands  of  variables  and  constraints.   However,  the  spareness 
of  the  tableau,  the  0-1  values  of  most  of  the  variables,  and  the  parallelism 
of  programming  algorithms  indicate  areas  for  future  development  of  the  inte- 
ger programming  approach  for  the  site  location  problem. 
3.2  A  Zero-One  Formulation 

It  is  well  known  that  any  arbitrary  integer  programming  formulation 
can  be  converted  to  a  zero-one  formulation  by  the  appropriate  addition  of 
zero-one  variables  for  each  nonzero-one  variable  [12].   The  number  of  new 
variables  added  can  be  controlled  if  knowledge  can  be  obtained  about  the 
upper  bounds  for  nonzero-one  variables  which  will  still  permit  all  required 
feasible  solutions.   The  purpose  for  reformulating  the  site  location  problem 
from  an  all-integer  to  a  zero-one  problem  is  to  examine  the  computational 
differences  between  the  two  major  integer  programming  solution  techniques. 
A  cutting  plane  algorithm  was  explored  for  the  all-integer  formulation.   In 
this  section  an  implicit  enumeration  algorithm  for  solving  zero-one  formula- 
tions will  be  examined. 


3.2.1  Formulation  of  the  Initial  Tableau 

Only  a  few  changes  have  to  be  made  to  equation  sets  (l)  through  (7) 

of  section  3.1-1  in  order  to  convert  to  the  zero-one  framework.   First, 

equation  set  (k)   can  be  omitted  since  the  zero-one  formulation  implies  that 

all  the  activity  variables  must  be  zero-one.  All  that  remains  is  to  convert 

the  utility  activity  variables  U. 's  into  expressions  of  zero-one  variables. 

To  illustrate  how  the  transformation  is  made,  suppose  U.  has  an  upper  bound 

of  B.   Then  whenever  U.  appears  in  the  formulation  it  is  substituted  by  the 

equivalent  binary  representation 

U.=  1  U.  +  2  U.  +  k   U.  +...+  2h_1  U. 
Xl  X2     X3  Xh 

where  U.   is  zero-one  and  h  is  such  that  2  -1  >  B. 

1 .  — 

J 

It  is  important  to  find  a  small  upper  bound  B  so  that  the  number 
of  variables  added  to  the  formulation  is  minimized.  This  is  accomplished  by 
returning  to  the  section  3.1.1  to  study  the  arguments  which  established 
equation  set  (3)» 

U.-  U.+  nL. .  <  n-1,  for  all  L. .  in  the  graph. 

It  was  shown  that  an  upper  bound  on  the  U. 's  can  be  determined  by  numbering 
the  nodes  of  a  solution  tree.  A  smaller  upper  bound  can  also  be  determined 
by  the  same  arguments.   Consider  any  solution  tree.   The  root  of  the  tree 
is  defined  to  be  position  1  and  each  node  is  given  the  position  number  ac- 
cording to  its  height  or  level  from  the  root.   Figure  15  displays  this 
numbering  scheme.  Many  nodes  will  have  the  same  position  number  which  was 
not  the  case  in  section  3.1.1.  All  nodes  not  in  the  solution  tree  are  given 
the  position  number  zero.   By  following  the  exact  arguments  of  equation  set 
(3)  it  can  be  seen  that  the  maximum  position  number  which  represents  the 
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Figure  15.   Another  Node  Numbering  Scheme  for  Solution  Trees. 
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maximum  solution  tree  height  is  an  upper  bound  for  each  U. .  This  is  true 
whenever  the  tree  height  is  less  than  n,  the  total  number  of  nodes  in  the 
region,  which  is  certainly  a  safe  assumption. 

It  is  important  to  note  that  the  maximum  tree  height  depends  on 
the  size  (number  of  parcels)  of  the  site,  the  size  of  the  geographic  region, 
and  the  particular  parcel  configuration  of  the  region.   For  example,  in  a 
6x6   rectangular  grid  parcel  configuration  the  maximum  tree  height  B  for  a 
site  of  maximum  size  of  p  parcels  is  given  by 

B  =  minimum  ( | p/2 | ,  12 ) 
where   f    means  the  integer  not  less  than  f .   This  would  be  an  area  site 
whose  parcels  took  the  shape  of  a  coiled  path  with  the  root  of  the  tree  at 
the  center  of  the  site.  For  most  practical  problems  you  can  reduce  B.   In- 
deed, restriction  of  the  tree  height  to  be  less  than  the  maximum  required 
could  be  used  as  a  device  to  control  such  pathologically  shaped  optimal  sites 

Table  k   shows  the  canonical  zero-one  formulation  of  the  site 
location  problem.   Formulas  for  the  number  of  variables  and  constraints  re- 
quired to  formulate  a  zero-one  area  location  problem  in  an  p  x  p  rectangular 

parcel  configuration  are: 

2 
number  of  variables  =  (5+h)p  -  ^p 


and 


2 
number  of  constraints  =  7P  -  6p  +  k . 


The  sparceness  is  estimated  by 

(28+8h)p2-  (2l++8h)p        .,,  -2, 
sparceness  =  * r— c =  0(hp   ). 

((5+h)p  -  Up)(7p  -  6p+U) 
where  h  is  the  number  of  binary  digits  required  for  each  utility  variable 
U.  to  be  expressed  in  binary  form. 
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OBJECTIVE  FUNCTION: 


n 

MINIMIZE 

Z 
i=l 

w. 

1 

N. 

1 

SUBJECT  TO  THE  CONSTRAINTS: 


FOR  DIRECTED  TREES 


(1)    Z  N.-      Z      L.   <  1; 
i=l    all  i,j  s.t. 
L.  .  in  graph 


Z       L.  .-   Z  N.  <_  -  1 
all  ij  s.t.   1J  i=l  X  " 
L. .  in  graph 


(2)       Z      L.k  <  1,  k=l,2:,...,n 

all  i  s.t. 

L.,  in  graph 
lk 


(3)    (1U.  +  2U.  +...+  2h"1U.  )  -  (1U.  +  2U.  +...+  2h'1U.  )  +  nL. .  <  n-1 , 
Xl    X2  Xh       h  J2  Jh      1J  " 

for  all  L. .  in  the  graph. 


(U)    2L. .+  2L..-  N.-  N.  <  0,  for  all  L..s  L..   pairs  in  the  graph. 


See  Table  2  for  additional  constraints  for  forming  an  area,  path,  or 
point  tableau. 


Table  k.      A  Canonical  Zero-One  Formulation  of  the  Site  Location  Problem. 
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3.2.2  Computational  Experience 

Computational  tests  were  run  on  small  uniform  grid  configurations 
analogous  to  the  testing  procedure  for  the  all-integer  formulation.  A  pub- 
lished FORTRAN  IV  code,  based  on  an  extension  of  the  implicit  enumeration 
methods  developed  by  Balas  and  Geofferion,  was  used  to  test  the  zero-one 
formulation  [l8].   Testing  was  conducted  on  an  IBM  System  360/75*   Weights 
were  arbitrarily  chosen  for  several  p  x  p  uniform  rectangular  grid  configura- 
tions. Figure  16  shows  the  average  results  of  several  problems  of  varying 
size.   Just  as  in  the  all-integer  formulation,  the  computational  time  is 
seen  to  be  directly  proportional  to  both  the  size  of  the  entire  region  as 
well  as  the  size  of  the  site  within  the  region.  Although  the  execution  of 
ALGOL  on  the  B67OO  and  FORTRAN  on  the  360/75  are  difficult  to  compare,  it 
can  be  seen  by  reviewing  figure  lU  that  the  execution  times  are  within  the 
same  order  of  magnitude. 

Several  schemes  were  provided  by  the  zero-one  code  for  assigning 
the  value  of  a  free  variable.  A  detailed  discussion  of  these  schemes  can  be 
found  in  [l8].   Each  scheme  had  a  profound  effect  on  the  execution  times  — 
sometimes  by  two  orders  of  magnitude  —  with  "scheme  two"  consistently  pro- 
viding the  best  results.   Schemes  that  take  advantage  of  the  site  location 
problem's  specific  structure  may  provide  future  efficient  algorithms. 

As  indicated  by  the  graph  in  figure  16,  the  computational  times 
explode  as  the  size  of  the  problem  expands.   Results  of  tests  on  a  6  x  6 
configuration  are  given  in  table  5  to  show  this  condition.   Note  that  the 
location  of  an  optimal  area  of  size  5  took  182  seconds  while  an  optimal 
area  of  size  10  was  not  found  even  after  75^+2  seconds  had  elapsed. 

It  is  interesting  to  observe  that  implicit  enumeration  algorithms 
will  determine  an  optimal  solution  before  it  can  indicate  it  is  optimal. 
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Figure  16.   The  Average  Processing  Times  of  Several  Zero-One  p  x  p  Problems 
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optimal  solution 
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was  found 
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iterations 

1 
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3 

1 

1 

73 
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12 

12 

10 

10 

311 
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15 
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13 

13 

383 
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6k 

6k 

^7 

kl 

1597 

5 

182 
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155 
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599i+ 

10 

75^2 

- 

- 

2U78 

300000 

15 

1895 

- 

- 

- 

75000 

20 

1789 

- 

- 

- 

75000 

Table  5.   Zero-One  Test  Results  on  a  6  x  6  Parcel  Configuration. 
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Also  suboptimal  feasible  solutions  can  be  indicated  as  the  algorithm  pro- 
ceeds.  Perhaps  these  early  results  can  be  explored  in  the  future  as  a 
heuristic  solution  to  the  site  location  problem. 
3.3  A  Zero-One  Multiple  Site  Formulation 

The  initial  tableau  for  the  multiple  site  problem  described  in 
section  2.k   is  constructed  by  concatenating  several  single  site  tableaus . 
The  interactions  between  each  site  are  expressed  within  the  objective 
function. 

Consider  the  following  integer  activity  variables  and  coefficients 

For  each  node  i,  1=1,2, ... ,n,  and  site  s,  s=l,2,...,m  let 

„s    fl,  if  node  i  is  used  for  site  s 
i    ) 0,  otherwise 

fl,  if  the  link  directed  from  node  i  to  node  j  for  site  s 
L. .  =  <    is  in  the  solution  tree  for  site  s 


\o,  otherwise 


h  -1 

US  =lUS+2U.+...+2S  U. 

\  X2  \ 

s 

=  the  binary  representation  of  the  utility  variable  for  node  i 

h 

used  for  site  s  where  h  is  such  that  2  s-l  >  B  ,  the  tree 

s  —  s 

height  for  site  s 

g 

U.   =  0,  1  as  determined  by  an  integer  programming  algorithm 
\ 

a.  =  area  of  parcel   i 

w?  =  weight (greater  than  zero)  of  using  parcel   i  for  site  s 

d. .  =  distance  from  parcel   i  to  parcel  j 

w  .  =  weight  of  using  parcel  i  for  site  s  and  parcel  j  for  site  t 


U8 
The  weight  w  .   could  be  determined  by  a  function  of  the  distance  d. .  and/or 

S"C1J  1J 

a  function  of  the  direction  with  which  parcel  j   is  situated  in  relation  to 
parcel  i . 

g 
A  =  lower  bound  on  the  area  for  site  s  if  site  s  is  an  area  site 

o 

A  =  upper  bound  on  the  area  for  site  s  if  site  s  is  an  area  site 
u 

A  canonical  formulation  for  the  multiple  site  location  problem  is  given  in 

table  6.   Note  that  the  objective  function  is  quadradic .  This  condition  can 

be  eliminated  by  the  following  well  known  transformation.   If  x,  y,  and  z 

are  zero-one  variables  then  z  =  xy  if  and  only  if  x+y  ^_  2z  and  x+y  <_  z+1. 

Therefore,  by  letting 


js  n1: 

stij   1  j 


Y__.  ,=  NT  N, 


where  Y  .  . .  is  a  zero-one  variable,  and  adding  the  two  constraints 
stij 

-is  -  nS  2r  . . .  <  o 

i    j     stij  — 
and  (6) 

NS  +  N*  -  Y  .  .  .  <  1 
l    j    stij  — 

to  the  tableau  of  table  6  for  each  Y  . . . -  the  multiple  site  problem  will  be 

stij 

converted  to  a  conventional  zero-one  integer  programming  problem. 

In  its  current  state,  the  formulation  permits  two  or  more  sites  to 
partially  overlap  or  occupy  the  same  parcels.  This  might  be  desirable  for 
certain  applications.   If  necessary,  it  can  be  avoided  by  the  addition  of 

the  following  set  of  constraints. 

m   s 

E  N.  <  1  ,  i=l,...,n  (T) 

.   l  — 
s=l 

Each  inequality  allows  a  specific  parcel  to  be  used  by  at  most  one  site. 
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OBJECTIVE  FUNCTION: 


SUBJECT  TO  THE  CONSTRAINTS: 


m 

n 

s 
w. 

1 

NS 

1 

m 

m 

n 

n 

NS 

l 

J 

MINIMIZE 

Z 
s=l 

Z 
s=l 

+ 

Z 
s=l 

Z 
t=l 

Z 
i=l 

Z 

st 

ij 

FOR 

DIRECTED  TREES 

n 

(1) 

Z   NS  -      Z       LS.  <_  1; 
i=l   X   all  i,j  s.t,   1J 
L?  .  in  graph 

s     n   s 
Z       L. .  -   Z  N.  <-l  ,  for  each  site  s=l,2,...,m 

all  ij  s.t.   J   i=l 

L?  .  in  graph 

(2) 

Z       L^k  <  1,  k=l,2,...,n,  s=l,2,...,m 

all  i  s.t. 

L?n  in  graph 
lk 

h  -1                   h  -1 

(3) 

(1US  +  ...  +  2  S  US   )  -  (1US  +  ...  +  2  S  US   )  +  nLS.  <  n-1, 

Xl                                \                  h                                Jh       1J  " 

s                       s 

for  all  L. .  in  the  graph,  s=l  ,2,  .  .  .  ,m 

(h) 



2  LS.+  2  LS.-  NS-  NS  <  0, 
ij      Ji    i    J  - 

s    s 
for  all  L..,  L..  pairs  in  the  graph,  s=l ,2, . . . ,m 

-1- J      J  -L 

Table  6.   A  Canonical  Zero-One  Formulation  of  the  Multiple  Site  Problem, 
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FOR  AN  AREA  FOR  SITE  s 


n 

(5a)    E  a.NS  <  AS  ; 
i  i  —  u  ' 
i=l 


n 

I      (-a.)NS  <  -AS 
i=l    X   X  ~   £ 


FOR  A 

PATH  FROM  WS 

l 

TO  NS  FOR  SITE  s 

(5b) 

all 

E 

k  s.t. 

<LL»  i  1  « 

lk 

in  graph 

all 

E 

k  s.t. 

(-Lik'  i  -1  » 

lk 

in  graph 

NS 
J 

li; 

-j 

1  -i 

FOR  A  POINT  FOR  SITE  s 


n 
(5c)   E  &*<.!; 

i=l 


n 

E   (-NS)  £  -1 
i=l    X 


Table  6  (continued) 
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Formulas  follow  for  the  number  of  variables  and  constraints  re- 
quired to  construct  a  linear  zero-one  multiple  site  location  in  a  p  x  p 
rectangular  parcel  configuration.   The  formulas  are  based  on  a  model  for 
multiple  nonoverlapping  area  sites.   This  includes  the  equation  sets  (l) 
through  (5a)  of  table  6  and  equation  sets  (6)  and  (7)  above. 

2  h        2   2 

number  of  variables  =  m((5+h)p  -  Up)  +  (p  -  P  )(m  -  m) 

=  (m  -  m)p  +  (6m  +  mh  -  m  )p  -  (^m)p 

2  h        2   2        2 

number  of  constraints  =  m(7p  -  6p  +  k )   +  2(p  -  p  )(m  -  m)  +  p 

=  (2m2-  2m)p  +  (9m  -  2m2+  l)p2  -  (6m)p  +  km 

where  m  is  the  number  of  sites  and  h=h  ,  8=1,2, ... ,m. 

The  formulas  do  not  include  those  variables  and  constraints  for  Y  ,  . .  where 

stij 

w  .  .  .=  0  when  s=t  or  i=j . 
stij 

It  is  interesting  to  note  that  if: 

(1)  only  point  sites  are  to  be  located  which  eliminates 
the  need  for  equation  sets  (l)  through  (U); 

(2)  the  number  of  points  to  be  located  is  n  which  is  equal 
to  the  number  parcels  in  the  region; 

(3)  no  overlap  is  permitted  (one  point  is  to  be  assigned 
to  each  parcel  in  the  region); 

(k)     the  integer  and  zero-one  requirements  for  variables  are 
dropped 
then  the  model  for  multiple  site  location  becomes  identical  to  the  quadratic 
assignment  problem: 

n   n  n   n   n   n 

minimize      Z        Z     w.  N.  +  Z        Z        Z        Z     w  ^ . .  N.  N. 
,  .  ,   11     .  .  .  .  _  .  n   stij   1   1 
s=l  1=1         s=l  t=l  1=1  j=l     ° 
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n 


g 

subject  to:    E  N.  =  1,  s=l,2,...,n 
i=l   X 


n 

Z  NS  =  1,  i=l,2,...,n 
s=l  x 


This  mathematical  programming  problem  was  first  formulated  by  Kocpmans  and 
Beckman  for  the  determination  of  optimal  plant  location  assignments  where 
costs  of  transportation  flows  were  to  be  considered  [19] •   It  is  also  being 
studied  by  Ray  as  part  of  an  urban  activity  assignment  model  [20].   The  ad- 
vantages of  the  model  presented  in  this  thesis  over  the  quadratic  assignment 
model  are  that : 

(1)  areas,  paths,  and  points  can  be  considered  simultaneously; 

(2)  high  resolution  planning  is  allowed  because  more  than  one 
parcel  can  form  a  site; 

(3)  structure  of  sites  are  very  general  due  to  the  linear  graph 
definition  of  sites; 

(h)     the  problem  can  be  reduced  from  quadratic  to  linear  form 
because  of  the  zero-one  formulation. 

Computational  tests  were  not  run  on  the  multiple  site  problem  due 
to  its  excessive  size  compared  with  the  single  site  formulations.   Indeed, 
it  is  apparent  from  the  results  on  single  site  problems  that  the  algorithms 
used  in  this  thesis  are  totally  impractical.   Suggestions  for  future  research 
into  finding  acceptable  solutions  to  the  site  location  problem  are  presented 
in  the  next  section. 
3.1+  Suggestions  for  Future  Research 

Four  possible  categories  of  future  research  into  the  development 
of  economic  site  location  algorithms  are: 
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(1)  zero-one  algorithms  taylored  to  the  specific  structure 
of  the  site  location  problem; 

(2)  heuristic  solutions  to  the  integer  programming  formulation; 

(3)  other  heuristic  techniques; 

(k)      fourth  generation  parallel  processors. 

Section  3*2.2  discussed  the  computational  effect  of  different 
schemes  for  assigning  values  to  free  variables  in  implicit  enumeration  al- 
gorithms.  The  suggestion  was  to  explore  schemes  that  take  advantage  of  the 
particular  formulation  presented  in  this  thesis. 

Because  of  the  extensive  size  of  the  problems  formulated  in  this 
thesis,  it  appears  that  suboptimal  approaches  to  the  integer  programming 
formulation  would  hold  the  most  promise  in  the  immediate  future.   At  the 
least  it  is  possible  to  accept  the  best  feasible  solution  produced  at  a 
specific  cut-off  time  by  an  implicit  enumeration  algorithm. 

Techniques  other  than  integer  programming  can  be  explored.   A 
simple  heuristic  algorithm  is  presented  in  the  next  chapter  for  area  loca- 
tions.  The  algorithm  treats  geographic  parcels  and  their  weights  as  a 
3-dimensional  topography  and  seeks  the  valleys  formed  by  low-weighted  par- 
cels as  possible  suboptimal  areas.   Critical  path  algorithms  could  be  used 
in  the  case  of  path  locations.   This  can  be  done  by  shifting  the  weights 
from  nodes  to  links  in  the  directed  linear  graph  model  of  a  parcel  system. 
The  weight  of  each  link  oriented  into  a  node  would  be  the  weight  of  that 
node. 

All  the  possibilities  that  have  been  explored  or  suggested  so  far 
have  an  inherent  parallel  nature  which  can  be  exploited  by  parallel  proces- 
sors such  as  the  ILLIAC  IV.   Chen  studied  an  ILLIAC  IV  version  of  the 
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implicit  enumeration  algorithm  used  in  section  3.2.2  [21].   It  was  shown 
that  ILLIAC  IV  can  solve  zero-one  problems  with  a  speed-up  of  Ik   times  that 
of  serial  computers  which  does  not  take  into  consideration  the  additional 
speed  of  the  ILLIAC  logic  circuits.  The  tests  for  the  zero-one  formulation 
were  run  on  a  360/75  which  can  execute  approximately  one  million  instructions 
per  second  (l-mip).   An  ILLIAC  IV  processing  element  has  been  estimated  to 
execute  3-mips.   Therefore,  we  could  expect  a  speed-up  of  approximately  k2 
times  when  solving  site  location  problems  on  ILLIAC  IV.   Chen  also  proposes 
a  special  purpose  processor  for  solving  zero-one  formulations  which  further 
reduces  computational  times. 
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CHAPTER  k 
LANGUAGE  CONSTRUCTS  FOR  EXPRESSING  PARCEL  WEIGHTS 

Most  users  of  geographic  based  information  are  not  computer 
scientists.  For  this  reason  generalized  geographic  information  systems  have 
been  developed  which  provide  users  with  simple  but  powerful  problem  oriented 
languages  in  order  to  retrieve,  analyze,  and  display  spatial  data.   Since 
these  users  possess  the  knowledge  required  to  meaningfully  weight  parcels  for 
site  location  problems,  it  is  important  to  explore  those  language  constructs 
which  allow  them  to  most  readily  create  and  use  parcel  weights. 

This  chapter  describes  the  language  extensions  made  to  the  geo- 
graphic information  system  NARIS  (Natural  Resource  Information  System)  to 
provide  non-programmers  with  a  capability  to  create  weighting  expressions. 
Weighting  expressions  in  their  simplest  form  are  analogous  to  PL/1  or  ALGOL 
arithmetic  expressions  except  that  identifiers  are  replaced  by  a  limited  set 
of  data  element  names  occurring  in  the  geographic  data  base. 

In  order  to  test  weighting  expression  constructs  for  locating 
optimal  sites,  a  small  interactive  system  was  developed  employing  a  sub- 
optimal  site  location  heuristic.   This  system  will  be  discussed  and  demon- 
strated with  actual  data  processed  by  NARIS  weighting  expressions. 
k.l     An  Introduction  to  NARIS 

NARIS  is  a  computerized  geographic  information  system  which  was 
developed  at  the  Center  for  Advanced  Computation  of  the  University  of 
Illinois  at  Urbana-Champaign  in  cooperation  with  both  the  Northeast  Illinois 
Natural  Resource  Service  Center  and  the  Northeastern  Illinois  Planning 
Commission  under  a  Ford  Foundation  grant  [8,22].   NARIS  is  designed  to 
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provide  inexpensive,  easy  to  use,  interactive  information  retrieval  capabil- 
ities to  agencies  and  individuals  in  need  of  geographically  locatable  data. 
A  brief  description  of  the  system  follows  in  terms  of  its  principle  data 
structures  and  the  major  conversational  language  components  provided  to  users. 
A  complete  syntax  of  the  language  can  be  found  in  Appendix  C. 

Currently,  each  parcel  in  the  system  geographically  corresponds 
to  the  legal  land  description  of  the  Public  Land  Survey  and  is  one  quarter- 
quarter  section  of  land  {UO   acres).   The  Public  Land  Survey  was  intended  to 
produce  a  uniform  grid  parcel  configuration.   However,  a  small  percentage 
of  parcels  are  oddly  shaped  or  missing  due  to  survey  errors  (c.f.  figure  l). 
The  resulting  configuration  must  necessarily  be  handled  by  generalized  parcel 
techniques. 

Figure  IT  displays  the  NARIS  data  base  structure.   The  primary 
units  of  data  are  called  data  elements  which  are  names  for  the  actual  data 
entries  in  the  data  base.   All  data  elements  are  aggregated  into  data  class- 
es.  In  the  figure,  "forestry"  is  a  data  class  with  data  elements  "type", 
"acres",  and  "density".   Data  element  values  may  either  be  numerical  quanti- 
ties or  character  codes.   In  the  forestry  example,  values  for  forestry  type 
are  the  character  codes  "pine",  "oak",  etc.   The  values  for  forestry  acres 
are  numerical  quantities  ranging  from  "00. l"  to  "UO.O".   In  the  event  that 
the  value  of  an  element  is  not  recorded,  a  special  "missing  data"  value  is 
stored. 

A  single  class  entry  in  the  data  base  is  called  an  occurrence.  An 
occurrence  has  one  set  of  data  values  for  each  data  element  within  the  data 
class.   For  example,  an  occurrence  of  forestry  data  would  be  the  data  element 
values  describing  a  single  homogeneous  stand  of  forestry.   The  stand  might 
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Figure  IT.   The  MARIS  Data  Structure 
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result  in  the  data  occurrence:   "type"  equals  "pine";  "acres"  equals  "15. 7"; 
and  "density"  equals  "0.75" •   Although  the  structure  for  each  occurrence  is 
fixed,  there  can  be  an  arbitrary  number  of  occurrences  for  each  class  in  a 
parcel. 

The  scope  of  a  retrieval  must  be  specified  before  a  request  can  be 
executed.   Scope  is  defined  in  NARIS  by  specifying  a  region.   Internally,  a 
region  is  simply  a  list  of  record  numbers  representing  a  subset  of  parcels. 
Regions  may  be  created  by  users  in  three  ways: 

(1)  by  explicit  geographic  description  using  standard  Public 
Land  Survey  notation  to  address  one  or  many  parcels; 

(2)  by  the  data  content  of  parcels  using  Boolean  expressions 
to  specify  the  conditions  that  data  values  must  satisfy; 

(3)  by  constructing  new  regions  from  existing  regions  using  the 
set  operations  "exclude",  "union",  and  "intersect". 

Regions  perform  a  dual  purpose.   First  they  allow  a  user  to  manip- 
ulate data  using  a  geographic  notion  that  is  ingrained  in  their  thought 
processes.   Second,  they  provide  the  associative  keying  capability  necessary 
for  on-line  information  retrieval.   By  limiting  each  request  to  the  parcels 
specified  in  a  region,  serial  searching  of  the  entire  data  base  is  eliminated. 

In  order  to  calculate  weights  for  each  parcel  in  a  region,  weight- 
ing or  arithmetic  expressions  have  been  designed  to  mirror  the  NARIS  data 
structure.  At  the  lowest  level,  weighting  expressions  provide  the  ability 
to  arithmetically  combine  values  of  numerical  data  elements  within  the  same 
occurrence.   By  combining  low  level  arithmetic  expressions  with  Boolean  ex- 
pressions, a  single  value  for  a  class  can  be  computed  on  a  subset  of  the 
occurrences  in  a  parcel.  A  single  value  or  weight  for  a  parcel  can  then  be 
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computed  by  arithmetically  combining  the  values  calculated  for  different 
classes . 

A  retrieval  request  specifies  a  region,  the  specific  primarys  to 
be  retrieved,  and  the  action  to  be  executed  upon  retrieval  primaries.   A 
retrieval  primary  can  either  be  a  single  data  element,  an  arithmetic  ex- 
pression of  data  elements,  or  a  list  of  these  primaries.   The  actions  each 
have  a  unique  verb  in  the  request  language:   TABULATE,  CALCULATE,  MAP,  and 
OUTPUT.   TABULATE  lists  information  on  a  parcel-by-parcel  basis.   CALCULATE 
displays  summaries  of  data  over  the  complete  region.   MAP  produces  scaled 
and  accurate  maps  of  parcels  that  are  gray  level  shaded  according  to  the 
computed  value  of  an  arithmetic  expression.   OUTPUT  produces  a  disk  file  con- 
taining the  value  of  one  or  more  arithmetic  expressions  for  each  parcel. 
The  OUTPUT  verb  allows  other  programs  such  as  report  generators,  statistical 
systems,  simulators,  and  site  location  algorithms  to  manipulate  the  data 
outside  of  NARIS  in  ways  that  would  be  beyond  the  sophistication  of  most 
users. 

Another  facility  in  the  language,  an  ABBREVIATION,  permits  a  user 
to  abbreviate  portions  of  a  request  with  a  single  name.   This  gives  the 
language  a  text  replacement  macro  capability  which  permits  long  requests  to 
be  broken  into  smaller  segments  and  also  lets  nonsophisticated  data  users 
easily'  execute  requests  created  by  technical  specialists.   English  text 
descriptions  of  data  classes,  elements,  and  values  are  available  through 
the  WHAT  IS  request.   This  mechanism  also  allows  users  to  recall  the  defini- 
tion of  an  entity.   Other  facilities  in  the  language  are  used  to  manage  the 
various  regions,  abbreviations,  and  other  entities.   For  example,  these 
facilities  allow  a  user  to  save  an  entity,  destroy  an  entity,  or  control 
which  users  and  agencies  have  access  to  an  entity. 
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4.2  Out putting  of  Arithmetic  Expressions  in  NARIS 

The  current  language  constructs  in  NARIS  for  specifying  weighting 
expressions  and  outputting  of  disk  files  of  parcel  weights  are  presented  in 
terms  of  their  syntax  and  semantics.   Future  language  concepts  will  be 
discussed  in  section  4.4. 
4.2.1  Fundamental  Definitions 
4.2.1.1  Syntax 

<empty>:  :=  a  null  string 

<letter>:  :=  A  |B  |.  .  .  |Z 

<digit>:  :=  0  |l  . .  .  |9 

<special  characters  :=   <single  space>|/|' 

<character>:  :=   <letter >  |  <digit  >  |  <special  character;- 1. 

<sign>::=  +|-|<empty> 

identifier >:  :=   <letter>  |<identifier><letter>  |<identifierxdigit  > 

<character  code>::=   <character>  |  <character  codexcharacter> 

<unsigned  integer >: :=   <digit> | <unsigned  integer ><digit> 

<integer>::=  <sign><unsigned  integer> 

<unsigned  number >: :=  <unsigned  integer> | . <unsigned  integer> 
|<unsigned  integer>.<unsigned  integer > 

<number>::=  <sign><unsigned  number> 

<class  identifier>: :=  <identifier> 

<data  element  identifier>:  :=  <identifier>  |  <class  identifierxidentifier> 

<numeric  data  element  identifier>: :=  <data  element  identifier 

<data  element  value>: :=  <integer> | <number> | <character  code> | "<character  code> 
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4.2.1.2  Semant  ic  s 


A  <class  identifier >  is  the  name  of  a  data  class  in  the  data  base. 
A  <data  element  identifier>  is  the  name  of  a  data  element  of  one  of  the 
classes.   A  <data  element  value>  is  one  of  several  possible  values  for  a 
particular  <data  element  identifiers   A  <data  element  identifier>  can  have 
either  numerical  or  character  values  but  not  both.  A  <numeric  data  element 
identifier>  is  a  data  element  which  has  <integer>  or  <number>  values.   Each 
reference  to  a  data  identifier  in  a  NARIS  request  implies  the  necessary  re- 
trieval of  the  corresponding  data  values  from  the  data  base. 

Some  of  the  <data  element  values>  will  contain  <special  characters> 
These  values  must  be  enclosed  in  quotes  (")  in  a  NARIS  request. 
4.2.1.3  Examples 

Consider  table  7  which  shows  three  possible  occurrences  of  soils 
data. 


SOIL 

NUMBER 

ACRES 

SLOPE 

W103 

15 

A 

27 

20 

C 

W103 

5 

B 

Table  7«   Some  Occurrences  of  Soil  Data 


SOIL  is  a  <class  identifiers   SOIL  NUMBER,  ACRES,  SLOPE  are  examples  of 
<data  element  identif ier>s .   One  possible  <data  element  value>  for  SOIL 
NUMBER  is  W103.   SOIL  ACRES  is  a  <numeric  data  element  identif ier>  because 
all  of  its  values  must  be  <number>s. 
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U.2.2  Functional  Transformations 

At  times,  the  arithmetic  operations  of  addition,  subtraction, 
multiplication,  division,  and  exponentiation  will  not  be  adequate  for  per- 
forming a  desired  calculation.  For  this  reason  functions  are  permitted  in 
arithmetic  expressions. 
I+.2.2.1   Syntax 
<function  request>::=  <the  part>  FUNCTION  <function  identifier><is  part> 

<function  definition^ 
<function  identifier>: :=  <identifier> 
<f unction  definition>: :=  <coordinate  pair><coordinate  pair> 

|<function  definition><coordinate  pair> 
<coordinate  pair>::=   (<coordinate><limit  indicator> ,<coordinate> ) 
<coordinate>: :=  <number> 
< limit  indicators  :=  +|-|<empty> 
<the  part>::=  THE | <empty> 
<is  part>::=   IS|<empty> 
it. 2. 2. 2   Semantics 

A  function  request  is  used  to  define  a  function  to  NARIS.   The 
function  is  used  within  an  arithmetic  expression  by  specifying  its'  identify- 
ing name.   The  coordinate  pairs  of  a  function  definition  represent  the  end- 
points  of  the  line  segments  of  a  piecewise  linear  function  of  one  variable. 
The  defining  coordinate  pairs  must  be  listed  by  giving  the  argument  (x-coord- 
inate)  first  and  the  function  value  (y-coordinate=f (x) )  second.   Coordinate 
pairs  must  be  listed  so  that  a  function  is  unambiguously  defined  by  specify- 
ing the  x-coordinates  in  increasing  order  from  left  to  right. 

The  limit  indicator  is  used  to  specify  discontinuous  functions.  The 
"+"  sign  after  the  coordinate  "x"  in  the  coordinate  pair  (x+,y)  states  that 
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the  function  value  "y"  is  not  given  to  the  argument  value  "x"  hut  to  the 
argument  value  just  larger  than  "x".   Similarly,  the  x-coordinate  "x-" 
indicates  a  number  just  less  than  "x". 
U.2.2.3  Examples 

Figure  18  shows  the  graph  of  a  function  given  the  identifier 
TREECUTCOST .   TREECUTCOST  represents  the  cost  of  removing  forestry  from  a 
parcel  of  land.   The  statement  for  defining  TREECUTCOST  is 

FUNCTION  TREECUTCOST  IS  ( 0,0 )( 50 ,100) ( 50+, 150) (100,200 )# 
U.2.3  Arithmetic  Expressions  Within  an  Occurrence 
U.2.3.1  Syntax 
<occurrence  arithmetic  expressions :=   <occurrence  term> 

| <unary  operator><occurrence  term> 

|<occurrence  arithmetic  expression><add  operator><occurrence  term> 
<unary  operator>::=  +|- 
<add  operator>::=  +|- 
<occurrence  term>::=  <occurrence  factor> 

|<occurrence  termxmultiply  operatorxoccurrence  factor> 
<multiply  operators  :=  *|/ 
<occurrence  factor>::=  <occurrence  primary> 

|<occurrence  factor><exponentiate  operatorxoccurrence  primary> 
<exponentiate  operator>::=  ** 
<occurrence  primary>::=  <unsigned  number> |<occurrence  variahle> 

|<occurrence  function> | (<occurrence  arithmetic  expression>) 
<occurrence  variable>::=  <numberic  data  element  identifier> 
<occurrence  function>::= 

FUNCTI0N[<function  definition> ] (<occurrence  arithmetic  expression>) 

|<function  identifier>(<occurrence  arithmetic  expression>) 


a 


50  100 

PERCENT    OF    FORESTRY 


Figure  18.   The  Graph  of  the  Function  TREECUTCOST, 
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U.2.3.2   Semantics 

An  occurrence  arithmetic  expression  is  a  calculation  rule  to  be 
executed  on  the  numerical  values  of  an  occurrence  primary.  The  value  of 
a  number  is  itself.   The  value  of  a  numeric  data  element  identifier  is  the 
data  element  value  represented  by  the  identifier  of  the  data  occurrence 
under  evaluation.   For  occurrence  functions,  the  value  is  the  result  re- 
turned by  the  transformation  defined  by  the  function.   Finally,  the  value  of 
an  occurrence  arithmetic  expression  enclosed  in  parentheses  is  the  result  of 
the  arithmetic  calculation  defined  by  the  recursive  analysis  of  values  of 
primaries  of  the  other  three  types.   All  integers  are  converted  to  reals 
and  all  arithmetic  is  real.   The  precedence  of  arithmetic  operations  are 
defined  by  the  syntax. 

It  is  important  to  note  that  the  expression  is  evaluated  only  on 
the  data  within  the  same  occurrence  which  is  data  that  necessarily  belongs 
to  the  same  class.   An  occurrence  arithmetic  expression  can  not  be  composed 
solely  of  numbers.  At  least  one  numeric  data  element  identifier  must  be 
included  within  the  expression  and  at  least  one  data  element  identifier 
must  include  a  class  identifier. 
U.2.3.3  Examples 

Consider  the  class  FORESTRY  with  numeric  data  elements  ACRES, 
TOTALPARCELACRES,  and  DENSITY.   Some  examples  of  occurrence  arithmetic  ex- 
pressions are: 

(1)  FORESTRY  ACRES 

(2)  (  FORESTRYACRES  /  TOTALPARCELACRES  )  *  100.00 

( 3 )  TREECUTCOST ( (FORESTRYACRES/TOTALPARCELACRES )*100 ) 

where  TREECUTCOST  is  the  function 
identifier  defined  in  section  U.2.2.3 
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(U)   -DENSITY**FUNCTION[ (0,0) (10. 0,0) (20. 0,5)] (DENSITY) 
4.2.U  Boolean  Expressions  Within  an  Occurrence 

Each  parcel  -will  probably  have  several  occurrences  of  data  for 
any  particular  class.  Arithmetic  expressions  can  be  defined  that  perform 
calculations  on  the  values  of  an  occurrence  arithmetic  expression  evaluated 
on  a  portion  of  the  data  occurrences  within  a  parcel.   In  order  to  specify 
which  occurrences  should  be  included  in  a  higher  calculation,  a  Boolean 
facility  is  provided  which  returns  a  TRUE,  FALSE,  or  MAYBE  value  for  each 
occurrence  in  the  parcel.  MAYBE  values  occur  when  relations  can  not  be  re- 
solved due  to  missing  data  or  execution  time  errors  such  as  divide  by  zero. 
U.2.U.1   Syntax 
<occurrence  Boolean  expressions :=  <occurrence  Boolean  term> 

|<occurrence  Boolean  expression>0R<occurrence  Boolean  term> 
<occurrence  Boolean  term>::=  <occurrence  Boolean  factor> 

|<occurrence  Boolean  term>  AND  <occurrence  Boolean  factor> 
<occurrence  Boolean  factor>::=  <occurrence  relation> 

| (<occurrence  Boolean  expression>) 
<occurrence  relation>::=  <data  element  identifier><equality  relational> 
<data  element  value> 

| <data  element  identifier><one  of  relational>(<data  element  value  list>) 

|<occurrence  arithmetic  express ion ><relational> 
<occurrence  arithmetic  expression> 

|<occurrence  arithmetic  expression><relational><number> 

| <number><relational><occurrence  arithmetic  expression> 

|<occurrence  arithmetic  expressionxone  of  relational>(<integer  list>) 
<equality  relational^  :=  IS|lS  N0T|<is  part>  EQUAL  <to  part> 
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|<is  part>  NOT  EQUAL  <to  part>|DOES  NOT  EQUAL 

| EQUALS | EQ | NE | NQL | NEQ 
<is  part>::=  IS|<empty> 
<to  part>::=  TO|<empty> 

<one  of  relational>:  :=   <is  part>  ONE  OF|<is  part>  NOT  ONE  OF 
<relational>:  :=   <equality  relational> 

|<is  part>NOT  <gtr  lss  partxthan  part> 

|  <is  partxgtr  lss  partxthan  part> 

|  <relational>  OR  <relational> 
<gtr  lss  part>::=  GREATER | GT | GTR | LESS | LT | LSS 
<than  part  > : : =  THAN | <empty> 
<data  element  value  list>::=   <data  element  value> 

|<data  element  value>,<data  element  value  list> 
<integer  list>::=   <integer> | <integer>,<integer  list> 
U .2 . U .2   Semantics 

An  occurrence  Boolean  expression  is  a  rule  for  computing  a  logical 
value  over  the  data  within  one  occurrence.  The  rules  of  evaluation  are 
analogous  to  the  rules  for  evaluating  arithmetic  expressions.   Occurrence 
relations  are  given  the  value  true  if  the  relation  is  satisfied  for  the 
arithmetic  expressions  or  data  elements  involved.   If  the  relation  is  not 
satisfied  then  the  value  is  false.   If  the  relation  is  indeterminable  because 
of  missing  data  or  execution  errors  then  the  value  is  maybe.   The  precedence 
of  Boolean  operations  are  defined  by  the  syntax.   The  meaning  of  the  logical 
operators  "AND"  and  "OR"  are  defined  in  Table  8. 

Comparing  the  value  of  an  arithmetic  expression  to  an  integer  list 
is  accomplished  by  rounding  the  value  to  the  "entier (value  +  0.5)"- 
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OPERANDS 

OPERATIONS 

bl 

b2 

bl  AND  b2 

bl  OR  b2 

true 

true 

true 

true 

true 

false 

false 

true 

true 

maybe 

maybe 

true 

false 

true 

false 

true 

false 

false 

false 

false 

false 

maybe 

false 

maybe 

maybe 

true 

maybe 

true 

maybe 

false 

false 

maybe 

maybe 

maybe 

maybe 

L 

maybe 

Table  8.   The  Meaning  of  the  Logical  Operators  AND  and  OR 
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1+.2.U.3  Examples 

Consider  the  class  FORESTRY  with  numerical  data  elements  ACRES 
and  DENSITY  and  the  character  data  element  TYPE.   Some  examples  of  occurrence 
Boolean  expressions  are: 

(1)  FORESTRY  TYPE  IS  PINE 

(2)  FORESTRY  ACRES  IS  GREATER  THAN  OR  EQUAL  TO  15 

(3)  FORESTRY  TYPE  IS  ONE  OF  (OAK,  MAPLE,  PIKE) 

(k)      (FORESTRY  ACRES  *  DENSITY)  LT  10.0  AND  (TYPE  IS  NOT  PINE  OR 
ACRES  GT  5) 
i+ . 2 .  5  Arithmetic  Expressions  Within  a  Class 
U.2.5.1  Syntax 
<class  arithmetic  expression> : := 

<functional  designator><occurrence  arithmetic  expression> 
|<functional  designator><occurrence  arithmetic  expression> 
<occurrence  where  clause> 
<occurrence  where  clause>::=  WHERE  <occurrence  Boolean  expression> 
<functional  designators  :  = 
SUM | TOTAL 
| AVE | AVERAGE 
| MIN | MINIMUM 
| MAX | MAXIMUM 
| COUNT 
U  .2  . 5 .2  Semantics 

A  class  arithmetic  expression  is  a  rule  for  computing  a  single 
result  on  the  data  within  some  of  the  occurrences  in  a  class  in  one  parcel. 
A  vector  of  values  is  produced  by  the  evaluation  of  the  occurrence  arithme- 
tic expression  on  each  qualifying  data  occurrence  within  each  parcel  record. 


TO 
If  a  where  clause  is  not  expressed  then  every  occurrence  within  the  class  is 
a  qualifying  occurrence.   When  a  where  clause  is  included  in  the  class  ex- 
pression the  Boolean  expression  is  evaluated  on  each  occurrence.   If  result- 
ing value  is  false  then  the  data  occurrence  does  not  qualify,  if  the  value 
is  true  then  the  occurrence  qualifies,  and  if  the  value  is  maybe  then  the 
computation  is  terminated  on  the  tract  and  the  user  is  notified. 

The  resulting  vector  of  values  are  operated  upon  by  the  algorithm 
associated  with  the  specified  functional  designator.   SUM  or  TOTAL  causes 
the  final  result  of  the  class  arithmetic  expression  to  be  the  summation  of 
values  in  the  vector.   AVERAGE  or  AVE  computes  the  average  value.   MINIMUM 
or  MIN  and  MAXIMUM  or  MAX  finds  the  minimum  or  maximum  value  respectively. 
COUNT  returns  the  number  of  values  in  the  vector  (equal  to  the  number  of 
qualifying  occurrences). 

If  a  class  expression  contains  a  where  clause  which  when  evaluated 
on  a  parcel  shows  that  none  of  the  occurrences  qualify,  then  a  value  of  zero 
is  assumed  for  the  result  regardless  of  the  functional  designator.   This 
condition  will  be  noted  by  the  system  and  will  be  indicated  to  the  user  if 
necessary. 
U.2.5.3  Examples 

Again  it  should  be  noted  that  all  the  data  elements  specified  in 
a  class  arithmetic  expression  must  belong  to  the  same  data  class.   Consider 
the  class  POPULATION  with  numeric  data  elements  SIZE  and  AVESALARY  and 
character  data  element  RACE.   Some  examples  of  class  arithmetic  expressions 
are: 

(1)  TOTAL  POPULATION  SIZE 

(2)  MIN  (POPULATION  SIZE  *  AVESALARY) 
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(3)   TOTAL  POPULATION  SIZE  WHERE  RACE  IS  SPANISHAMERICAN 
(U)   MAXIMUM  (POPULATION  AVESALARY  *  (SIZE  **  2)) 

WHERE  (SIZE* AVESALARY)  LE  10.0  AND  RACE  IS  NOT  CAUCAISAN 
h.2.6     Arithmetic  Expressions  Over  a  Parcel 
U.2.6.1   Syntax 

<parcel  arithmetic  expressions :=   <parcel  term> 
| <unary  operator><parcel  term> 

|<parcel  arithmetic  express ion >< add  operator><parcel  term> 
<parcel  term>::=  <parcel  factor> 

|<parcel  term><multiply  operator><parcel  factor> 
<parcel  factor>::=   <parcel  primary> 

|<parcel  fact or >< exponentiate  operator><parcel  primary> 
<parcel  primary>::=   <unsigned  number > | <par eel  variable> 
|<parcel  function> | (<parcel  arithmetic  expression>) 
<parcel  variable>::=   <class  arithmetic  expression> 
<parcel  function>::= 

FUNCTION  [<function  definition>]  (<parcel  arithmetic  expression5, ) 
|<function  identifier>(<parcel  arithmetic  expression>) 
U.2.6.2   Semantics 

A  parcel  arithmetic  expression  is  a  rule  for  computing  a  single 
result,  which  could  represent  a  weight,  based  on  all  the  data  within  a 
single  parcel.   The  rule  expressed  by  a  parcel  arithmetic  expression  is 
executed  on  the  numerical  values  of  parcel  primaries.   The  rules  are  com- 
pletely analogous  to  the  rules  of  occurrence  arithmetic  expression.   The 
value  of  a  number  is  itself.   The  value  of  a  parcel  variable  is  the  result 
of  the  class  arithmetic  expression  calculation  on  the  data  element  values 
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of  the  specified  class.  For  parcel  functions,  the  value  is  the  result  re- 
turned by  the  transformation  defined  by  the  function.  Finally,  the  value 
of  a  parcel  arithmetic  expression  enclosed  in  parentheses  is  the  result  of 
the  arithmetic  calculation  defined  by  the  recursive  analysis  of  values  of 
primaries  of  the  other  three  types. 

If  a  class  arithmetic  expression  primary  of  the  parcel  expression 
does  not  contain  a  class  identifier  then  the  class  is  assumed  to  be  the 
class  of  the  class  arithmetic  expression  preceding  it.   In  order  to  determine 
the  class  of  the  first  class  arithmetic  expression  used  in  the  parcel  ex- 
pression, the  first  data  element  identifier  in  the  class  expression  must 
include  a  class  identifier. 
k.2.6.3     Examples 

Consider  again  the  example  classes  and  data  elements  defined  earli- 
er.  Examples  of  parcel  arithmetic  expressions  follow: 

(1)  10**2+5.0 

(2)  TOTAL  FORESTRY  ACRES 

(3)  (TOTAL  FORESTRY  ACRES )/( TOTAL  SOIL  ACRES )*100 
(h)      (MAX (POPULATION  SIZE*AVESALARY)WHERE 

RACE  IS  NOT  CAUCAISAN)  -  MIN  AVESALARY 
(5)   (EC0L0GYWT( TOTAL  FORESTRY  ACRES) )/(T0TAL(EC0N0MICC0ST( FORESTRY 

ACRES ))WHERE  DENSITY  GTR  5-0) 

where  EC0L0GYWT  and  EC0N0MICC0ST  are  previously 

defined  function  identifiers 
U.2.7  The  Output  Request 
U.2.7.1  Syntax 
<output  request>::=   <for  clause>  OUTPUT  ON  FILE  "<file  name>" 
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<parcel  arithmetic  expression  list># 
<for  clause>::=  FOR  <region  expression> 
<region  expression>: :=   see  Appendix  C 
<file  name>::=   <identifier> 
<parcel  arithmetic  expression  list>::=   <parcel  arithmetic  expression> 

|<parcel  arithmetic  expression  list> ,<parcel  arithmetic  expression> 
U .2.7 .2   Semantics 

The  output  request  provides  a  user  the  capability  of  forming  a  file 
of  calculations  on  NARIS  data.   A  record  is  produced  for  each  parcel  in  the 
region.   Each  record  contains  a  geographic  locator  and  the  values  of  each 
arithmetic  expression  written  in  "E"  format.   The  file  can  be  read  by  pro- 
grams other  than  NARIS. 
U.2.7 .3  Examples 

An  example  is: 

FOR  COOKCOUNTY  OUTPUT  ON  FILE  "WEIGHTS" 
(ECOLOGYWT( TOTAL  FORESTRY  ACRES) )/ (TOTAL 
(ECONOMICCOST( FORESTRY  ACRES))  WHERE 
DENSITY  GTR  5-0),  TOTAL  FORESTRY  ACRES# 

where  EC0L0GYWT  and  EC0N0MICC0ST  are  previously  defined 
function  identifiers. 
k.3     A  Weighting  Expression  and  Heuristic  Area  Location  Example 

In  order  to  explore  the  use  of  weighting  functions  for  site  loca- 
tion problems  over  regions  with  large  numbers  of  parcels,  a  suboptimal  heur- 
istic area  location  algorithm  was  developed. 
U.3.1  A  Heuristic  Area  Location  Algorithm 

The  central  idea  behind  the  algorithm  is  to  view  a  parcel  system 
and  its  parcel  weights  as  a  three-dimensional  topography.  The  x,y-dimensions 
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represent  the  geographic  coordinates  while  the  z-dimension  represents  the 
elevation  which  is  established  by  the  weights  of  each  parcel.  The  "valleys" 
of  this  topography  are  incrementally  "flooded"  in  order  to  locate  suboptimal 
low-weighted  areas.  At  each  level  of  flooding,  every  "puddle"  is  examined 
in  order  to  determine  if  any  satisfies  all  the  specified  constraints.  The 
most  important  constraints  will  be  those  that  dictate  the  final  surface  area 
of  a  puddle  which  has  the  same  size  as  the  site  represented  by  the  puddle. 
If  no  puddle  satisfies  all  the  constraints  then  the  flood  level  or  watertable 
is  increased  enlarging  each  puddle  and  merging  many  together.   The  algorithm 
terminates  when  one  or  more  puddles  at  the  present  watertable  satisfies  all 
the  constraints.  The  puddle  whose  weights  sum  to  the  minimum  and  satisfies 
all  the  constraints  represents  the  suboptimal  area. 

Adjacencies  are  determined  by  a  table  that  represents  the  linear 
graph  model  of  the  parcel  configuration.  Puddles  are  created  by  referring 
to  this  table.   Since  a  linear  graph  model  is  incorporated  into  the  algorithm, 
adjacencies  can  be  specified  or  ignored  without  regard  to  geographic  parcel 
contiguity.  Figures  19  and  20  shows  the  flowcharts  of  the  heuristic  area 
location  algorithm  and  the  puddle  location  subalgorithm.  Appendix  B  contains 
code  for  the  interactive  heuristic  program. 
1+.3.2  A  Housing  Subdivision  Example 

The  example  problem  is  to  suboptimally  locate  a  site  for  a  housing 
subdivision.  The  region  under  study  is  Marengo  township  located  in  McHenry 
County  within  the  State  of  Illinois.   This  township  is  composed  of  576  par- 
cels of  land  in  a  rectangular  grid  configuration  where  each  parcel  is 
approximately  one  quarter  mile  square  (U0  acres).  The  only  constraint 
will  be  on  the  number  of  parcels  that  compose  the  site. 
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Figure  19.   The  Flowchart  of  the  Heuristic  Area  Location  Algorithm, 
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Figure   20.      The  Flowchart   of  the  Puddle  Location   Subalgorithm 


77 

Actual  data  from  the  NARIS  data  "base  will  "be  used  to  establish 
each  parcel's  weight.   Detailed  descriptions  of  the  data  element  values  used 
in  the  weighting  expression  will  not  be  given.   This  information  can  be 
obtained  directly  from  NARIS.   Only  a  general  description  of  each  part  of 
the  weighting  expression  is  necessary  to  understand  what  the  request  is  do- 
ing. 

The  abbreviation  text  replacement  capability  in  the  NARIS  language 
will  be  used  to  segment  the  weighting  expression  in  order  to  facilitate  the 
explanation.   An  abbreviation  request  allows  a  user  to  specify  a  single  word 
that  represents  an  arbitrarily  long  text  string.   The  word  can  be  used  in  any 
request  by  following  it  with  a  period  ( . ) .   The  system  expands  the  text  as 
it  compiles  the  request.   For  example, 

THE  ABBREVIATION  ABBR  IS  ABBREVIATION  # 
will  allow  a  user  to  specify  "ABBR."  instead  of  having  to  type  out  the  word 
ABBREVIATION. 

The  object  of  the  weighting  expression  is  to  establish  a  weight  for 
each  parcel  which  represents  the  parcel's  relative  suitability  to  be  included 
as  part  of  a  housing  subdivision.   Each  segment  of  the  expression  will  weight 
a  parcel  high  if  it  is  very  suitable  for  housing  with  respect  to  the  data 
under  consideration  and  lower  if  it  is  not  suitable.   Since  the  heuristic 
site  location  algorithm  finds  minimally  weighted  sites,  the  final  weights 
will  be  multiplied  by  -1  in  order  to  make  the  least  weighted  parcels  more 
desirable. 

The  following  function  will  be  used  repeatedly  in  order  to  negativ- 
ly  weight  parcels  with  poor  characteristics. 

THE  FUNCTION  IGNOR  IS  (0,0) (0+,-1000) (100,-1000 )# 


78 
A  parcel  will  receive  an  extremely  low  weight  whenever  this  function  en- 
counters an  argument  whose  value  is  greater  than  zero. 

The  following  abbreviations  will  make  up  the  final  weighting 
function. 

THE  ABBR.  FLOODWT  IS  (IGN0R( TOTAL  STREAMS  FLOODACRES) )# 
Any  parcel  that  is  subject  to  flooding  will  be  ignored. 

THE  ABBR.  CONSTRUCTWT  IS  (IGN0R( TOTAL  GEOCONSTRUCT 

ACRES  WHERE  TYPE  IS  NOT  ONE  OF  (G2,  G3)))# 
Any  parcel  which  contains  poor  geological  conditions  for  construction  must 
be  ignored. 

THE  ABBR.  BADLANDUSEWT  IS  ( IGNOR ( TOTAL  LANDUSE1 

ACRES  WHERE  CODE  IS  ONE  OF  (33.0,  6k. 0,   11.0,   72.0,  9^.0) ))# 
Mines,  cemeteries,  sewage  disposal  works,  refuse  dumps,  and  marshlands  are 
extremely  undesirable  and  any  parcels  containing  them  will  be  ignored. 

THE  ABBR.  SOILWT  IS  (2*  TOTAL  SOIL  ACRES  WHERE 

SLIGHTSEPTICLIMIT .  AND  SLIGHTURBANLIMIT . )# 
Parcels  are  weighted  porportional  to  the  amount  of  soil  with  only  slight 
urban  and  septic  limitations.   SLIGHTSEPTICLIMIT.  and  SLIGHTURBANLIMIT.  are 
abbreviations  created  by  a  U.S.  Department  of  Agriculture  Soil  Conservation 
Service  scientist. 

THE  ABBR.  OPENLANDWT  IS  (2*  LANDUSE  ACRES  WHERE  CODE  IS  8l.0)# 
Open  land  is  desirable . 

THE  ABBR.  FORESTRYWT  IS  (80*  FUNCTION  [  (0,0)  (10,1 )  (1+0,-2)  ] 

(TOTAL  FORESTRY  ACRES ) )# 
Parcels  with  10.0  acres  of  forestry  are  best.   Parcels  with  less  may  require 
tree  planting  for  aesthics  and  those  with  more  may  require  clearing. 
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THE  ABBR.  DRAINAGEWT  IS  (2*  TOTAL  GEOWATER  ACRES 

WHERE  TYPE  IS  Gl)# 
Good  drainage  characteristics  are  desirable.   The  final  expression  is  defined 
by: 

THE  ABBR.  SUBDIVISIONWTS  IS  ■ 

-(FLOODWT.  +  CONSTRUCTWT.  +  BADLANDUSEWT .  +  SOILWT.  + 

OPENLANDWT.  +  FORESTRYWT .  +  DRAINAGEWT.)  +  273# 
■where  the  factor  273  is  an  adjustment  so  that  the  lowest  or  best  parcel 
weight  is  zero.   The  final  weights  are  output  to  a  disk  file  by 

FOR  MARENGOTWP  OUTPUT  ON  FILE  "SUBDIVISION"  SUBDIVISIONWTS.# 
The  constraints  on  the  number  of  parcels  for  the  site  are  constructed  from 
a  file  in  which  the  coefficient  "l"  is  output  with  each  parcel's  geographic 
locator. 

FOR  MARENGOTWP  OUTPUT  ON  FILE  "PARCELS"  1  # 

After  the  weights  have  been  output,  the  heuristic  program  is  in- 
voked.  Figure  21  shows  the  conversation  between  the  user  and  the  algorithm. 
A  colon  is  used  to  signal  the  user  that  a  response  is  expected.   A  map  of 
each  parcel's  weight  is  output  to  the  station  as  a  2k   x  2h   matrix  of  three 
digits.   Three  asterisks  (***)  represent  a  digit  greater  than  999-   It  is 
interesting  to  note  the  asterisk  pattern  outlines  a  stream  and  its  flood 
banks.'  The  next  map  shows  all  the  puddles  found  when  a  solution  was  detected. 
The  final  maps  display  all  the  suboptimal  areas. 
k.h     Future  Language  Concepts 

The  weighting  expressions  provided  in  NARIS  form  a  preliminary 
approach  to  the  problem  of  providing  a  comprehensive  user  oriented  language 
for  generating  parcel  weights.   Additional  language  and  system  concepts 
which  require  future  research  are  briefly  discussed  in  this  section. 


80 


RUN  AREA/HEURISTIC/OBJECT; END 
17:25??////  STAT10N=SELECTRIC  OF  AREA/HEURIST1C/0BJECT=2987  OfENED  //// 


2987  ARF.A/HEURISTIC/OBJECT@50*BOJ 
WHAT  IS  YOUR  NARIS  USER  NUMBER? 
:U004 


WHAT  IS  YOL'K  FILE  NAME  FOR  THE  PARCEL  WEIGHTS? 
: SUBDIVISION 


CONSTRAINTS  HAVE  THE  FOLLOWING  SYNTAX: 

?FILE  NAME??BLANKS??LE  OR  GE?? BLANKS?? NUMBER  WITH  A  DECIMAL  POINT? 


WHAT  IS  YOUR  FIRST  CONSTRAINT? 
: PARCELS  GE  20.0 


WHAT  IS  YOUR  NEXT  CONSTRAINT?   TYPE  NO  IF  YOU  HAVE  NO  MORE  TO  ENTER. 
:PARCELS  LE  24.0 


WHAT  IS  YOUR  NEXT  CONSTRAINT?   TYPE  NO  IF  YOU  HAVE  NO  MORE  TO  ENTER. 
:NO 


THE  PARCEL  WEIGHTS  ARE: 


185 193 190 196 194 179 189****** 1931 31*************** 1621 36 139 138 152 1281 70189 
187***193195195193*********175******113113103105131  93142155161109151140 
******  179 1642 10*********191 102******1 131 13  87145193159195207176237167132 
*** 193348207*********195 103******1 151 13 1151401 75 185167249200202 136226273 
*** 193 196********* 195 105  77***115115134176186155235227195268269249273*** 
*** 193198  lb  1******143 115****** 11 3153193 1831 7 1165 1421 77 193259273199409265 
200185143*********113111***173149157197185175195194113139182236217353209 
206  193*********115122******202202191195195176236155185203265235191273244 
*********191  99115122***233228197195195196195165123273265261211231433321 
*** 179 193135 115  95******176170194169193163156233191*********236206131385 
i72185149119  7  3  71***133170193237***175192369273******234171199273425433 
206119  51  46  54******129181191211***172146189241******181289195257329194 
94  69  55  57******  74132193154171***154156149******140155245233183188224 
65  50  81******  48  64  83137117129***117113************158186180123***210 
62*********  57  37  33  43  52*********  33*********  33***  73127163194***198 
*********  42  48  44  46  54  48  42***  58***************  44  42  52  72  66  82  78 
************  85***  0*************************************************** 
************************************************************************ 
***************************  55  53  51***  35  37  47  £,2  37******113  35****** 
************************  43  33  39  37  39  45  51  g2  73***  39******  74  35  60 
173253112  69  65101  54143  72  47  70  85  50  50  46  55164***123143151******  93 
*********  93136113  85  85  71  82  57  57  54  54  51133188191179153113********* 
143***196175107  83105  85  45113113  49  39  33  95172173189200235241193149131 
195193193175111  81  67  43  57  49  39  47  79  69163174165149211273273223207154 


ALL  THE  AREAS  ARE: 


Figure  21.   The  Conversation  Between  a  User  and  the  Hueristic  Algorithm, 
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113113103105  .   93  .   .   .  109  .   . 

102  .   .  113113  87  

103  .   .  115113115  

105  77  .  115115  

115  ..  113  

113111 113 

115 

.   .   .   .   99115  

.   .   .   .  115  95  

.   .   .  119  73  71  

.  119  51  46  54  

94  69  55  57  .   .   7A  

65  50  81  .   .   48  64  83  .  117  .   .  117113  

62  .   .   .   57  37  33  43  52  .   .   .   33  .   .   .   33  .   73  

.   .   .   42  48  44  46  54  48  42  .   58 44  42  52  72  66  82  78 

.   .   .   .   85  .    0  


55  53  51  .   35  37  47  62  37  .   .  113  35  .   . 

43  33  39  37  39  45  51  82  73  .   39  .   .   74  35  60 

112  69  65101  54  .   72  47  70  85  50  50  46  55 93 

.   93  .  113  85  85  71  82  57  57  54  54  51 113  .   .   . 

.   .  107  83105  85  45113113  49  39  33  95  

.   .  Ill  81  67  43  57  49  39  47  79  69  


AREAS  THAT  MEET  CONSTRAINTS  ARE: 


115 

.   .   .   .   99115 
.   .   .   .  115  95 
.   .   .  119  73  71 
.  119  51  46  54  . 
94  69  55  57  .   . 
65  50  81  .   .   . 
62 


SUM  OF  WEIGHTS  =       1605 


TYPE  YES  IF  YOU  WANT  TO  KNOW  WHICH 

CONSTRAINTS  ARE  VIOLATED  FOR  AREAS  NOT  IN  THE  SOLUTION. 

TYPE  NO  OTHERWISE. 

:NO 


PROCESSOR  TIME  =     7  SECS. 


Figure  21  (continued) 
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i+.i+.l  Distance  and  Direction 

Calculating  a  distance  from  a  parcel  is  tedious  in  the  current 
NARIS  system.  A  data  class  called  MAPDATA  with  data  elements  X  and  Y  is 
part  of  the  data  base.   Each  occurrence  of  MAPDATA  is  an  X  and  Y  coordinate 
of  a  corner  point  of  the  polygon  defining  the  parcel.  The  center  of  a  parcel 
is  approximated  by  the  average  X  and  Y  coordinates  of  the  parcel.  The  follow- 
ing expression  uses  the  formula  for  calculating  distance  between  two  coordin- 
ates.  Suppose  the  distance  from  coordinate  (123,1+56)  is  to  be  formulated. 
The  expression  for  accomplishing  this  is 

((((AVE  MAPDATA  X)-123)**2)+( ( (AVE  MAPDATA  Y)-h^6)**2) )**0 . 5 
It  is  apparent  that  this  is  not  user  oriented.  Distance  should  be  a  primi- 
tive construct  within  any  generalized  geographic  information  system.  Within 
the  NARIS  syntax,  it  could  be  treated  as  an  addition  to  <parcel  primary>s 
defined  in  section  4.2.6.1.   The  syntax  could  be 
<parcel  primary>::=  <unsigned  number >  |  <par eel  variable2* 

|<parcel  function> 

|(<parcel  arithmetic  expression>) 

|<distance  specif ication> 
<distance  specifications  :  = 

<the  part>  DISTANCE  <from  part><parcel  specif ication> 
<from  part>: : =  FROM | <empty> 
<parcel  specif ication>: :=   <x,y  coordinate> 

|<Public  Land  Survey  parcel  specification> 
<x,y  coordinate >: :=   <number>,<number>  |  <number  number> 
The  distance  expression  in  the  example  would  become 
THE  DISTANCE  FROM  123,^56 
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The  above  construct  could  easily  be  added  to  the  NAEIS  language. 
However,  distance  from  a  collection  of  parcels,  such  as  a  NARIS  region, 
presents  a  significant  problem.   The  problem  arises  when  the  parcels  are 
not  all  contiguous.   In  this  case  the  definition  of  distance  can  be  formu- 
lated using  many  schemes.   Each  algorithm  must  be  studied  with  respect  to 
speed  and  user  acceptance. 

The  notion  of  direction  can  also  play  in  important  role  in  parcel 
weighting.   For  example,  being  down  wind  from  a  factory  would  have  an  effect 
on  the  optimal  location  of  a  housing  subdivision.  Direction  and  distance 
could  be  combined  together  to  form  a  polar  coordinate  factor  for  weighting. 
U . U . 2  Numerical  Associations  for  Non-numerical  Data 

Only  numerical  data  elements  can  be  manipulated  with  arithmetic 
operators.   In  order  to  increase  the  power  of  weighting  expressions,  a 
mechanism  for  attaching  numerical  values  to  the  character  values  of  a  non- 
numerical  data  element  should  be  provided.  A  possible  construct  could  take 
the  following  form. 
Association  request>::  = 

ASSOCIATE<association  identified  WITH  <data  element  identified 
AS  Association  list># 
Association  identifiers  :=  <identifier> 
Association  list>::=   <association> 

|  Association  list>  ,<association>  |  Association  list><association> 
Associations  :=   <data  element  value>:  <number> 
An  example  with  the  data  element  FORESTRY  TYPE  might  be 

ASSOCIATE  FORESTRYRANK  WITH  FORESTRY  TYPE 

AS  PINE:-1,  OAK: 5-5,  MAPLE: 2# 


Qh 
An  association  could  then  be  permitted  in  any  expression  in  which  a  character 
or  numeric  data  element  identifier  was  allowed.   For  example,  the  following 
could  he  legitimate  parcel  arithmetic  expressions. 

(1)  (TOTAL  FORESTRY  ACRES)* (AVERAGE  FORESTRYRANK) 

(2)  TOTAL  FORESTRY  ACRES  WHERE  FORESTRYRANK  GTR  0 
Associations  could  he  handled  within  the  NARIS  system  in  the  same  way  as 

the  other  entities  —  Regions,  Abbreviations,  and  Functions  —  are  processed. 
All  entities  have  a  structure  characterized  by  some  information  corresponding 
to  a  user  specified  identifier  which  is  analygous  to  the  structure  of  asso- 
ciations . 
U . U . 3  The  Influence  of  One  Parcel  on  Another 

Activities  and  conditions  that  occur  in  parcels  within  the  sur- 
rounding region  of  a  particular  parcel  may  have  an  effect  on  the  final  weight 
that  a  parcel  is  given.   The  ability  to  calculate  a  factor  for  a  parcel  based 
on  data  in  parcels  that  are  within  a  specified  distance  should  be  studied. 
One  approach  to  this  problem  would  be  to  allow  a  <parcel  primary>  of  the 
following  form. 
<parcel  primary>::=  <within  specificationxclass  arithmetic  expression> 

|etc.  see  section  k.2.6.1 
<within  specif ication>: :=  WITHIN  <unsigned  number>  MILES 
Some  examples  are: 

(1)  WITHIN  2  MILES  AVERAGE  FORESTRY  ACRES 

(2)  WITHIN  3.5  MILES  TOTAL  LANDUSE  ACRES  WHERE  TYPE  IS  FACTORY 
Another  important  concept  is  related  to  the  requirement  for 

"searching  upstream"  for  various  conditions.   The  purchaser  of  a  parcel  of 
land  might  wish  to  know  if  a  factory  is  polluting  the  stream  that  runs 
through  the  parcel.   This  desire  leads  to  the  notion  of  a  linked  data  base 
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where  data  elements  within  one  record  can  point  to  data  in  another.   Mechan- 
isms for  tracing  pointers  and  retrieving  information  along  the  way  could 
provide  a  powerful  search  mechanism. 
1+.1+.4  Optimal  Regions 

Section  U.l  provides  an  introduction  to  the  NARIS  concept  of  a 
region  .  Recall  that  a  region  is  a  name  associated  with  a  list  of  parcel 
numbers  which  can  be  created  by  geographic  location,  data  content,  or  set 
manipulations.   Since  an  optimal  site  is  also  defined  as  a  special  set  of 
parcels  it  would  be  natural  to  incorporate  site  location  within  a  NARIS 
region  request.   This  is,  of  course,  a  hypothetical  discussion  until  effec- 
tive site  location  algorithms  can  be  developed.   One  possible  form  of  an 
optimal  region  request  could  be 
<region  request>::=   <the  part>  REGION  <region  identifier> 

<is  partxregion  expression># 
<region  expressions  :=   see  Appendix  C 

|<location  type>  WHERE  <parcel  arithmetic  expression>  IS 

optimization  typexconstraint  list> 
<location  type>::= 

PATH  FROM  <parcel  specification>  TO  <parcel  specif ication> 

|AREA  BETWEEN  <unsigned  number>  AND  <unsigned  numberxarea  unit> 
<area  unit  > : : =  ACRES | SQMILES 
<optimization  type>: :=  MAXIMUM | MINIMUM 
constraint  list>::=  SUBJECT  TO  <constraint> 

| constraint  list>,<constraint>  |  <empty> 
<constraint>:  :=  <parcel  arithmetic  expressionxrelationalxnumber> 
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For  example,  the  housing  subdivision  problem  of  section  U.3  with  an  addition- 
al population  constraint  would  be  expressed  as 

THE  REGION  HOUSING  IS  AREA  BETWEEN  .800  and  960  ACRES 
WHERE  SUBDIVISIONWTS.  IS  MINIMUM  SUBJECT  TO  TOTAL 
POPULATION  SIZE  LESS  THAN  200  # 
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CHAPTER  5 
CONCLUSIONS 

The  linear  graph  integer  programming  formulation  of  the  site 
location  problem  provides  several  flexibilities  that  might  not  occur  with 
other  approaches.   First  of  all,  it  was  shown  that  path  and  point  location 
can  he  viewed  as  special  cases  of  the  more  general  area  location  problem. 
Another  important  characteristic  is  that  the  algorithms  are  system  independ- 
ent.  This  means  that  these  techniques  can  potentially  be  added  to  existing 
generalized  geographic  based  information  systems.   The  condition  that  sup- 
ports system  independence  is  the  fact  that  the  model  and  algorithms  are 
linked  to  the  geographic  structure  of  the  data  instead  of  the  particular 
geographic  coordinate  schemes  used  in  data  collection  and  referencing.   Only 
characteristics  of  adjacency  are  important. 

It  is  important  to  note  that  adjacency  can  be  arbitrarily  defined. 
Optimal  areas  were  defined  in  terms  of  a  linear  graph  which  was  not  required 
to  be  planar.  Two  parcels  that  are  not  necessarily  contiguous  can  be  linked 
together  for  special  applications.  For  instance,  in  urban  regions  it  may  be 
more  important  to  link  parcels  with  respect  to  transportation  accessability 
rather  than  just  adjacency.  Other  linkages  may  occur  due  to  communication 
capabilities. 

Just  as  important  as  adding  links  is  the  ability  to  break  links. 
For  example,  in  regions  of  highly  varying  topographies  two  adjacent  parcels 
could  be  separated  by  several  hundred  feet  of  average  elevation  which  would 
prohibit  the  two  parcels  from  being  part  of  the  same  site. 

All  of  the  examples  thus  far  have  dealt  with  two-dimensional  parcel 
spaces.   Ths  is  not  a  restriction  since  the  linear  graph  model  can  easily 
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incorporate  a  three-dimensional  volume  system  as  well  as  an  n-dimensional 
system  of  space  and  time  adjacencies. 

The  links  of  the  graph  in  the  integer  programming  model  have  asso- 
ciated variables  which  can  be  weighted  in  a  manner  similar  to  node  weighting, 
These  weights  could  represent  degrees  of  connectivity.  By  incorporating  the 
link  weights  into  the  objective  function,  site  locations  could  optimized  not 
only  by  parcel  weights  but  by  the  interactions  between  adjacent  parcels. 

It  was  also  shown  that  the  addition  of  other  constraints  is  a 
simple  matter  in  the  integer  programming  model.  This  feature  also  allows 
particular  parcels  to  be  forced  into  the  solution.  Their  respective  node 
variables  merely  have  to  be  set  equal  to  one  in  the  initial  tableau.  The 
same  is  possible  for  link  variables . 

It  was  demonstrated  that  the  shapes  of  optimal  areas  can  be  con- 
trolled. For  instance,  long  narrow  areas  would  be  quite  undesirable  for 
many  applications.   This  condition  can  be  eliminated  by  easily  implemented 
restrictions  on  the  height  of  the  solution  tree. 

A  major  feature  of  the  multiple  site  formulation  is  that  the 
seemingly  inescapable  quadratic  nature  of  the  problem  was  overcome.  The 
zero-one  variables  of  the  formulation  provide  a  key  to  reducing  the  problem 
to  a  conventional  linear  form. 

The  most  serious  problem  with  the  formulation  is  the  large  compu- 
tational times  required  to  solve  integer  programming  formulations.  Future 
research  is  required  into  optimal  and  suboptimal  techniques  that  capitalize 
on  the  unique  structure  of  the  models  presented  in  this  thesis.   Of  course, 
the  special  properties  of  paths  and  points  can  be  exploited  to  develop  lo- 
cation algorithms  more  efficient  than  integer  programming.   For  example, 
critical  path  techniques  could  be  used  to  locate  optimal  paths.   Location 
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of  an  optimum  weighted  point  would  be  trivial  since  a  point  is  defined  as  a 

single  parcel.   Only  a  simple  sort  has  to  "be  performed  on  a  list  of  parcel 
numbers  and  their  associated  weights.   Positional  relationships  are  not 
important  since  questions  of  contiguity  do  not  arise  when  searching  for 
points. 

User  oriented  language  constructs  for  expressing  parcel  weights 
were  demonstrated.  They  have  been  successfully  implemented  with  a  system 
that  has  gained  user  acceptance.   However,  several  concepts  are  still  lack- 
ing and  require  further  research. 
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aijRRnuc-HS    R6T00    ALGOL    COm»KEb    LEVEL    I  1 , 0? .  00, 0  1  5 ,       THURSDAY,    Ob/29/72,    1?iM    AM, 


OPTTMAl/aREA 


REr.lN 

itjiHU'uutunumiHni>mrMtu»tnm«ju»u»m»uutntntni 

»  TUNING  PARAMETERS 

MtinoitstituuiutinnuituiusinitituiiiitiitiiMtiiiMsntnti 

nCTTNE 

uUMIFNOOES" 
10 

fl 

t  nuhqE"  of  notes  tl,r,   parcels! 

nUM^FLInks* 

,t 

t. 

t    NUMBER  OF  LINKS 
..UmTtEn. 
* 

1. 

I    NUMBER  nr  ITERATIONS. .fACH  Kill  HAvE  a  otfeerfnt 
I  AREA  CONSTRAINT 
NUM"FCONSTRAlNT$. 
n 

f  , 

I  NUMBER  Of  constraints  OTHER  than  area  bounds 

uAX'ImEs 

i  eo">o 
«l 

1  MA„JMUM  CUMULATIVE  PPOCFSSHR  TImE  IN  SI»TYITmS  OF  A  SEC 
v/»l.'lE  ARRAY 

HlGHBOAREASt 
«•%. )•?»«. 6 

)  I 

*  UPPER  ROUNDS  PF  «OL.UTtON  AREAS  TNDEy  BY  ITERATION 
I  NUMBER 

>i*L'IE  ARRAY 

LO«RDARE*Sf 
a. 3. 1 .2,4,6 

)  I 

t  lOnEP  ROUNDS  HE  «OLnTtPN  AREA  InOEifEn 

*  BY  ITEPATlnN  NUM»ER 
Al  IE  ARRAY 

NO^EARf AS( 
5.2. i .2. 3,2,  1  .  1  ,  l  .2 
)  I 

I  AREA  Or  EACW  Nnnr  (DABCFL) 
vAL'iE  ARRAY 

«EK,hTs( 


oooiOOnoiO 

RLOCk  Is  SEGMENT  nOOn? 


00?:00n0 

0 

002tOOoO 

0 

00?:00n0 

0 

002i OOnOi 

0 

002  I'oonn 

0 

002i00n0 

0 

00210000 

0 

00?,00nn 

0 

00?| 00r>0 

0 

00?,00no 

0 

00?,00o0 

0 

00?, oono 

0 

00?,00on 

0 

00? 1OO00 

0 

00?  1OO0O 

0 

00?,00o0 

0 

00?,oono 

0 

00? lOOnO 

0 

00?iO0rn 

0 

o«2toeeo 

0 

oo?7oono 

0 

•02,0000 

0 

00?IO«00 

0 

oojioono 

0 

00?  ,00r>0 

0 

OOJtOOnO 

0 

oo?^ooon 

0 

06?i00n0 

0 

oo?,ooon 

10 

[)ATA  is  000* 

.0Nr 

oo?~oono 

0 

oo? joono 

0 

00?|0000 

10 

002ioonn 

10 

00?t00o0 

0 

0  0  ?  ,'  0  0  r\  n 

lO 

RATA  IS  000a 

.Our, 

oo?,oonn 

0 

00?,00n0 

|0 

0  0  ?  •'  0  C  n  0 

(1 

00?t00nn 

10 

0  0  ?  :  0  0  n  n 

10 

0  0  ?  •"  0  0  n  0 

lO 

nATA  IS  000» 

.ONr, 

00? .  ooon 

0 

00?,00o0 

|0 

00?,oono 

|0 

9U 


«.2.|.2.5,5.2.3.2.3 
)> 

»   «ek,ht  or  each  parcel 

IItfltt(*l(tt(tttlXt>I>ttttttfttf««>ftVltltttIt(TXlI1IllIttfItltt*lt«II>tI 

FlLF    LInECMYUSE^KInO'HS.A'JFFFPW.MAXrECsiZE^.INTMOOF'A  )l 

rur  c*Rn(KiN0.9.8uFrERSit2.HiXRrr^i?F.i« ,  i 

PEFTNE 

NuucoL^»(2*NUMnrNonf ^*NUMnrtlNK«»i >   «» 

I  NUMBERS  OF  ACTIVITY  VARIABLES  PLUS  J  FOR 
«  THE  RIGHT  HAND  S'OE  COLUMN 
Nu»R0"S»<5*2*NUN0FLINKS*NUM1FLINRS  fll  V  2»«»NUM0FN00ES 
♦NUMOrCONSTRAINTS)  »t 

*  NUMBER  OF  CONSTRAINT*  NEEOEp  TO  S*TISFT  ArEa 
t  AND  CONTIGUITY  RfSUIRFMrNTS 

INTEGER  ARRAY  A P Jf 1  I NUNOFNOOES» t INU"OFNOOES , , 

«  TELLS  IF  NODE  "FIRST  INf»E»"  IS  ADJACENT  TO 
«  NOoE  -SEcONn  TNBfl"  MMEN  ThF  V*l  lit  IS 
I    GREATER  THAN  JCRP.   THE 

*  VALUE  IS  THE  LIN'  NUMBER 
TArllNUMCOlS«0*NUMRO»Sl' 

nErtNE  At  I»jj«TAf j, i j  ,, 

I  THE  INTEGER  PR06R»MTNfi  TABLfAU  MATRf* 
INT'OER  I.JI 
fNTfSE"  TIM, 

t  PROCESSOR  TIME  IK  EACH  ITEBaTIOn 
INTEGER  ITERI 

«  ITERATION  NUMBER 
INTEGER  CUMTIMEl 

I  CUMULATIVE  PROCESSOR  TIME  Or  ALL  PROCESSING 

t  NOT  COUNTING  DATA  INITIALISATION 
INTEGER  lOTIMtl 

«  CUMULATIVE  I/O  TIME 
TNTfGER  PIvOTNUM, 

t  NUMBER  OF  INTEGER  PROGRAMMING  PTvOTS 
L*ifL  NOSOLUTION.STOPI 
LABFL  ne»ti 
nEFTNE  STEPl(I«L.U).FOR  ILL  STEP  1  UNTIL  U  00  ft 

nEFTNE 

STArTlINk*  NUMOrNOOF^»l  *• 

t  TABLEAU  COLUMN  KMERF  LINK  VARIABLES  START 
ENnLINM.NUHoFNnpES+NUMO'LlNKn  f, 

*  TABLEAU  COLUMN  WERE  LlN*  VARIABLES  END 
^TARTUT ILITY»ENPLTNK*1  '• 

t  TABLEAU  COLUMN  WHERE  UTILITY  VARIABLES  START 
ENnUTILlTV.NliMenLS»1  «l 

t    TABLEAU  COLUMN  W*ERF  uTTLITY  VARIABLES  ENO 
nFr'NE 

LlNK( I , j)«AOJC I , Jl  fl 

>  LIN«  NUHRfR  OF  LINK  RftlNTlNG  TO  NOOE  J  FROM  NODE  I 

nEFTNE 

STARTNP[)E»1  ft 

I  TABLEAU  COLUMN  WERE  NO*F  VARIABLES  START 
EN0N0DF»NU"0FNpnFs  *l 

I  TABLEAU  cnLUMN  WERE  NpnF  VARIABLES  ENO 


00, 
00? 
OAT*  IS 
002 
00? 
00? 
IS 
00? 
IS 


DATA 


DATA 


00? 
00? 

002 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00; 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 
00? 


lOOoO 

i  OOnO 

OOOi 

lOOnO 

1 0  0  n  o 

lOOnO 

0007 

i  0  0  n  n 

000* 

i'OOoO 

t  OOnO 

(ooon 
(oono 
t  00 no 
lOOoO 
t'OOftO 
iOOoO 
•oono 

jOOpJ 
,00n3 
ioons 
t'00p3 
lOOpJ 
j00n« 
i  oonft 
ioons 
ioonft 

lOOnft 

i  o  o  n  * 
iOOnft 
i00n6 
i00n6 
lOOnft 
lOOnft 
i00n6 
t00n6 
100ot 
|00p6 

iOOnt 
iOOnft 
t00p« 

,"  0  0  n  ft 
tOOpft 

•oonft 
1 0  0  n  6 
,'00n6 
lOOpft 
tOOnft 
)00n6 
,00rft 
lOOnft 
:O0D6 
lOOnft 
lOOnft 
lOOnft 
:00n« 
lOOnt 


o 

0 

Our, 

0 

0 

0 

ONf. 

0 

ONr, 

0 

0 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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nCFtNC 

itNKPOst j >«i*NuMorNnoES  *» 

t    TABLEAU    COLUMN    or    LIN„    T 
>tniiiri«itutunmitnmmiuiniuHtn«)tumiiini«tt>iUi>ti 

pRO'EDURE   CREATEADJI 

t  C"E4TES  TME  ADJACENCY  MATRIy  ADj  FROM  DAT*  CAPOS 
REGTN 
tNT'SER  ARRAY  Tf MPtOlB] I 

sTt*l(  It  ltNUNOr«JOOf  S) 
BEGIN 

READ(C»RD.«9!5>.STEBt(J»Oi85    TEHPfJl)> 
FUR    ji.l     STEP    i    UNTIL    7    00 
IF     TEMPtJl     OTR       0     'MEN 
A0J[TEHP[01,TEMPtjn»"TEMPtUlll 
ENP| 
fNO    CREATEADJI 

l«ItttltttIIttSISIIIIItttS«tttt««'tIX«t«IttltIISlltI(StIttIItItlIIISIISI 

PRO'EOUBE  FtLLNATRIXi 

t  THIS  IS  A  MOOFL  SPECIFIC  TABLEAU  MATPI*  GENERATOR. 

t  CONSTRAINTS  ARE  "ENFRATfD  In  CLASSES  THAT 

*  CORRESPOND  TO  THE  EOIIATTON  TYPES  SPECIFIED  FOR 

I  TNt  XODtL, 
RCITN 
tNT'GEB  R| 


PROCEDURE  lUMTxil 
I 

I  ONE  MORE  NODE  TH»N  LtNKl 
I 

t   N 

I  I—         t-- 
I   »  Nf J]  .   »   LfT.jj  .  , 
I  /•-        /•- 

I  I"!         ALL  E I • J  J  IN  GRAPh 
I 

rEgtn 

STE»1( J.STArTNODE.ENONOOE) 
BEGIN 

A(B. JJ l»"I I 
AC»»! . JJ l»l I 

ENni 

STEP1(  J»STA|»Tl.  iNBf  ENOLlNR  ) 
BEGIN 

A(Bi J) lal| 
ACP*l,j]i«»lt 

FNOI 
AtR.NUNf 0LSll»-1 I 

AIR*1 .NUHCDLS1 |a]| 
RI»B»2| 
FND  suBNTXli 


PROCEDURE  SUBMTK?i 
I 
X  DIRECTED  TBFFSl 


00?tOOnA 

00?iOOfiA 

0O?,00n*S 

00?,00n6 

00?i00n6 

oo?iOOn« 

oo?  toon* 

2 

O0?|00oft 

CREATFAOj  IS 

SEGmfNT  0 

00?i00r,2 

OOFiOOo? 

3 

OOriOOn? 

007,0015 

007,0016 

ooriooi  r 

oor,oom 

3 

oorioo?o 

CRE*Tt*DJ(007 

IS  00?A 

2 

00?|00(1» 

oo?Tooo« 

00?  iOOf.6 

00?, 0006 

00?i00n6 

00?|0004 

00?  lOOnft 

2 

00?i00n6 

FILLMATRIX  is 

SFGRfNT  0 

OOj,  tOOoC 

OOaiOOqO 

00a, 00(10 

oOjjOOrrO 

00a iOOoO 

00» lOOnc 

00a  i  OOnO 

OOtJOOno 

00»,00o0 

OOA^OOftO 

oo»teooo 

OOitOOpC 

3 

OOt Toono 

OOi^OOftO 

« 

ooiToooo 

oonoon2 

00»,00o5 

• 

oo4loonT 

OOi  lOOnr 

« 

00»  lOOnC 

oo»:oons 

OOAiOOnC 

* 

OOAiOOpF 

00.  .001  1 
0O»7O0l3 

00  j  ,-  00  1  « 

3 

00 » ,00) A 

00(  i  OOi a 

00 1  i  0  0 1  A 

0  0 1  i  0  0  i  • 

0 
0 
0 
0 
0 
0 
0 
0 

onr 

o 

o 

4 

2 

0 

0 

5 

1 

OnG 

0 

0 

c 

0 

0 

0 

0 

0 

004 

0 

0 

c 

0 
0 
0 
0 
0 
0 
0 
0 
0 

c 

0 

« 

5 

t 

3 

i 
2 

2 
5 

1 

0 
0 
3 

« 
« 

A 
4 
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i 

i  I  — 

I       >       LlI.Jl    LEO    1*    J«1.2. 

t    /"• 

t    *LL    Ltl»J]    Iw    GR'RH 

I 

«EGTn 
fNT'GER   tempi 


*TE»1< J.l.NuHOENOOtS) 
BEGIN 
5TEP1 U.I.NUHOFNODES5 

\f    TEHP|»L!NK(I.J)    «TR    0 

ThEN  »fR.LlNK'OS(TrMP)l»«ll 
»[R.NUmCOLS1I"1I 
RlaR»tl 
ENni 
fNO  <UBmTx2| 


PROCEDURE  SUBMTY3I 

t 

I  NO  DIRECTED  CIRCUIT*! 

I 

I  U(M"U(J]«N*L(I»-I]  LEO  N»l»  FOR  »LL  Ltl.jl  IN  GRAPH 

I 
*EfiTN 
tNTEGER  TEMPI 

sTE»l(I.l»NUMOrNOOES5 

STE'H J.l.NUMOENOOES) 

if     tempi«i.ink(  it  j^  otp  o 

THEN  BEGIN 

»tR.L!NKPOS(TrMP)l '.NUMprNOOESI 
»lR.ST»RTUTlLtTy*I-l )|«1 I 

»CR«ST»RTuTIlITy*J'1]I»"1| 
»tR.NUMCfLSl  l»NllMOENnOES»l  • 
R  l«R*l  I 
•     ENni 

FND  *URMTx3| 


PflO'EOURE  SUBNT*«| 
* 

I  NODE  VARJABELj,  »»E  Oil | 
» 

» ■  Nt  !  i    LEO  l,  I.I.? N 

* 
pEG'N 

<;TEe!  (J. STARTNODE -ENDNOTE) 
BEGIN 
»l  B. JJ  lal  I 
»t R,NUmCOLS]i«1i 

BlaRal I 

END  I 
FNO  SHBMTX6 1 


OOt  1 00 ) ft  I  4 
00»iOOi«i« 
OOl  i  00 l • I  * 

oo»;ooi»«« 

00tj00l4,4 
OOt iOOl»«* 

oo» Too i «i« 

3  00t"00l4i4 

SUBMTX2  IS  SEGMENT  OOOoB 
OORiOOQOlO 
OOBiOOOOlO 

«  OOglOOOOl* 

OOBiOOOOl* 
00p,,"0003l* 
00p,|0009l2 
OOPiOOoRlO 
00b,0OoCi2 

•  O0g,oopri* 

SUBMTX2(00B>  IS  QOli  long 
3  OO.tOOial* 

oo.,oo)4i« 

00»,0014|4 
00,700l4|4 
OOt |00t4l4 
00.  i'OO)  4  |4 
00»JOOl4|4 
00. J0014|4 

3  00>t00l4l4 

submtXs  is  sesmfnt  ooonc 

00c  lOOoO'O 
OOCiOOOOlO 
00r700nO|4 
00c  t'00o3l4 

4  00ci0004«3 
OOcjOOoTi I 
00C£000*|0 
OOciOOnO'O 

oociooor«5 

00C|00l ft, 1 
4  00Cl00l4t5 

SUBMTX3(00C)  IS  0017  LONG 
3  00»t00l4l4 

OOt  ,001 4  14 
00t,00|4|4 

00.(001*1* 

OOt t  0  0 1 4 «4 
OOt t 00 1 4 14 

OOt .00, «! 4 
00iJ00'l4|4 

3  0  0  i  ;'  0  0  1  4  >  4 
OOt (00)4l4 

4  00»i00l5]2 
OOt ,00171? 
OOt . 001 9  t  0 
OOt t  OOl 1 1  2 

a  OOt  :00l C I* 

3  00t:00in5 
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PROCEDURE  sum,,,  00*,00,C,5 

«  mk-lim.  rNrrp.rr.oN,  ••"••'"■'j 

*  l*Lri*Jl*f«LIJi|l-N|li*MIJI  LC«  O.FOR  All  LIt.Jl.LCJ.il  ool'oolcls 

*  PAIPS  IN  GRAPH  00»,00ici5 
aEGTN  00»|00ici' 

3  00a  1 00 1 C  15 

■  TE-|fI.l.NUllOMflDKll  SUBMTX5    ,S    "oSLn^00 

ir  km.hiiikcI.ji  ar,  „  SoS'SSSSlS 

»t«.LINKROS(L?HK(j.l, >,,.«,  SSn' SS2J  S 

«CR.STMTliaoE*I-t|l..||  DOnO     S 

»tR.ST»RTNOnr*J-lJl.-!l  OOn  00l1t« 

rini  0001001413 

CNO  SU8MTX5I  ^^^,"! 

*  ooniooi»'« 

SUBMTX5(000)  IS  001*  LONG 

RROtEOURC  SUBMTX6I  J              0^k'lVC,l 

„  00«,00iri3 

J  »»"  CONSTRAINT,  S'jJJJJjJ 

(   N  OO.tOOJc'S 

!  '"                 <""8  URRtRROUNO  "ijSSi^ 

«  >  (ARE»m.Ntin  <  ::  !::'!', 

«  '••               ■  ««  lowerppuno  ss;;ss;j;i 

i               "*  oo*7ootc«» 

BEGIN  OO^OOjClS 

sTEPlfJ.STAPTNOOE.ENONOOr)  3                                  22l!JI2!S!5 

PE6IN  OO*i00lC' 

A(R.JJ..N0ftC»RE»SfJ-STARTN0r>fli  ,                                2Sl!S«!«!5 

A{R»J.jJ,.-NO0E»REASfJ-«T*RTN0ntll  '                                  OoJjOOfS.J 

AtR.NU«cOLS)«"MI«Hgo»RfAStITt»-n«  »                                no!iS2!!lS 

.t.M.ju«eoiij...Lo.M«Kt.sifTc..,„  !!:;!!;;:? 

END   SUiMTK*,  oo.;oo»eii 

..........  oo»|00yn«* 

PROCEDURE    SUBHTxri  3                                  00Ai00»f>lS 

r  00»,00>ni» 

t    TO  makE  TA8LfAU  LEXtCBRAPHlr  PO.TTlVE.  Oojjoo'nl' 

t   N  00»j00?n«5 

n  1..  00»t00?ni5 

;   »,  ui«,  «.  «-.*«..,„/•  il:\ll$ll 

,  ,.,  oo4joo?ni5 

,  o0Ai00?nt5 

atGt^  00*T00>ril5 

sTE»lfJ.STARTUT!L!TY.ENOUTHTTY,  ,                                    n^1™'^ 

Afu.JII.li  3                                  00Ai00?m5 

00i,00?0i5 
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afRtNUMfCILSl  l"<*ilMOFN0OES*<N.IM0fP>'n0ES*l  ))    OTV    2    I 

Rl»»«l I 

rND  SUBmTx^i 

dRO'EOUrE  SUBMT*8I 

t 

»  TO  COMPLETE  TABLEAU  FpP  GOMnpri  i 

t 

I    nl  H  GEO  0.1.1.2 v 

t 

»  Ltl.jl  GEO  O.rnR  ALL  Lfl'Jl  lw  r.RAPH 

I 

tutll  GEO  0. t.l  .?•...  ,u 

t 
■  CQfN 

STE»1(J.1.NUMC0LS-I  ) 
BEGIN 

AtBi Jl la'll 
RlaB»l I 

Enoi 
f*d  SUBMTXBl 


l' 


pRO'EOURE  OBjECTi 
« 

*  OBJECTIVE  FUNeTfNi 
I 

I    M 
I   I--  |.. 

X      >      *E!GHTt  M«Nt  T  J    ♦      >      Ontlt.ll    ♦       >    O.Utll 

t    /••  /•*  /•• 

I    1»1  ALL    Lit. J)  T'l 

J  n     IRAPH 

» 

REGTN 

<TE»1 ( J.STARTNOPEtENONOnF) 

AtO.J}l»"EIGHTSrJ.START»OOEl| 

F*0    TRJFCTl 


pRO'EOUPE    ZfROAl 

1  CLEARS  OUT  TME  T«RLFAn 
=  EGtN 
^TE»](  !,O.M"iRO»S> 

STE°1( j.ltXUMCOLS) 

At I • Jl l "0 | 
rNO  ZE«nA| 

Plat  I 

7  F-  B  ">  A  i 
nBjFCTI 
tUR^T** I 
<  U  8  u  T  X  1  | 
CUBIT'S  I 
«;UB«T*4  I 
SU8UT  <2l 
<UB"TXF  I 
«:  U  B  ■•  T  X  3  I 

tuH-rxs  i 


00* 

i00l3l 

0 

00 

i  0  0 1 4 

5 

00. 

,0016 

1 

00/ 

t'0016 

2 

00/ 

|0016 

2 

00 

J0016 

2 

00/ 

100^6 

2 

OOi 

10036 

2 

00/ 

i'00*6 

2 

00 

,"0016 

2 

00 

,0016 

2 

00 

t'ooi* 

2 

00 

.00  16 

2 

00 

70016 

2 

00 

l'0016 

2 

00 

iToo^6 

2 

00 

U0016 

2 

00 

itOOlF 

0 

00 

(70019 

1 

00 

ll00l» 

3 

00 

LOOK 

} 

00 

uoom 

0 

00 

iiOOin 

0 

00 

>i'00in 

0 

00 

1  |00in 

0 

00 

1  joom 

0 

00 

iToo^n 

10 

00 

I  lOOlD 

0 

00 

A  t"  0  0  •*  n 

0 

00 

>tOO*.o 

0 

00 

k  t'OOin 

0 

00 

>  too^n 

0 

00 

uoo^n 

rO 

00 

lioom 

0 

0  0 

1  lOOin 

0 

00 

»,ooip 

0 

00 

1,00a? 

5 

00 

>i'00a3 

0 

00 

1  tOOo3 

0 

00 

1|00*3 

10 

00 

s  i00o3 

10 

00 

ll00«3 

0 

00 

I  .00a3 

|0 

00 

r,oo«3 

4 

00 

> ; OOaP 

0 

00 

l t'OOAR 

1 

00 

1  ,ooab 

1 

00 

■  , OOop 

s 

00 

1  iOOor 

3 

0  0 

1  ! 00«n 

1 

0  0 

I tOOoO 

!5 

00 

»  1 i  0  0  n  r 

.  3 

00 

>  :00ar 

1  1 

00 

itOOoF 

lb 

0  0 

I  [00*0 

1 3 

00 

1  tOOM 

11 

00 

Il00<1 

15 

9y 


FNO  FILLHATPJXI 

IlitttlKttSllStttXtlltttlltlttll'SttttttdtllltlttlllltMtltllltllttMt 

PROCEDURE  OUTPUTI 

»  OUT'UTS  THE  TABLEAU 
BEGIN 

DEFINE    H.NUMRCIHS    *•     N.NllMCOLS    (I 

INTEGER  ARRAY  1  NT  f 1  I6«] , 
INTEGER  I| 
FORMAT  Fl<X«.6»I2)l 
FORMAT  r2(I3'X] >6«T2) | 
f 

RRlTE(LINCfS"<I^  im 

FOR  li.I  STEP  1  UNTIL  N  0"  INt(!1|b!  OIV  1 0 1 

*RlTe(LlNE«fl 'FOR    T  I  •!     $TfP    t    UNTIL    N    00    InTMUI 

FOR    I'«t    STEP    1    UNTIL    n    po    INt;  I  )!•!•(( I    D1V    10>MO)1 

«RlTE(LINE,Ft .FOR  t|»l  STEP  1  UNTIL  N  00  INT[I])I 

WP.ITE(l!NE)I 

FOR  I»"0  STEP  J  until  M  pn 

KRITEfLINE.F?.  I.Fnp    Jl»l     STEP    I    UNTIL    N    Of)    ACI.JDI 
END    OUPUT, 

MtIttlI«tlS(tttlSI(I«tIttIStX«Itttttt«t«I*«t«II«II>ttltI«>IIt*IItttt«tt 

pROfEOURE  OUTPUTSOLNI 

t  OUTPUTS  NODES  IN  SOLUTION 
BEGIN 
INT'GER  TEMPI 

INTFGER  SUMI 

TEMPHNuMRONS-NUMCOLS  +  STABTNlOEl 
wRITE(LINE,<//. "SOLUTION  NOOTS  APF|"//»)I 
^TE»1( I. I iNuMOFNOOEs) 
BEGIN 

TEMPl«TtMR*l' 
JF  AtTE*P,NUNCOLSl  F0L  ' 
THEN  BEGIN 

wRITE(LI*E.<X10»"NrtOE  "»!?»«T)I 
SUMi«suM»NEIGMTSt I»n I 
ENPi 
END  i 
*RITE(LINE.</."SUM  OF  *EIGHT«»",T?>,SU")> 
FNO  OUTPUTSOLNI 

nnt*iiiuiiinninitnittiiiti«ntiuuttitnttuuiMMiiituuitii 

pRO'EDURE  GnMORTji 

t 

X     ALGORITHM  ?63»  MITH  WILSON'S  CUTS 

I  FROM, 

I  COLLECTEn  ALGOR'THM^  rRDM  CACM 

*  PAGES  263A-P1 .P?.P3 

I 
nEGTN 

DEFINE  M.NUMROWS  *•  N.N.IMCOLS  •• 

INTEGER  I.K.,i.L.R«C.T.S«LAMBfi»NI)N»  LAMppAOENOM  i 

BOOLEAN  NULL.HFLAGl 


00»i00^?l3 
EILLMATPI"(00A)  IS  OOfcP  LOnG 

2  00?t00n6l0 

00?tOOn*lO 
00?t00n6'0 
00?i00o6  10 

2  OOjiOOnAlO 

Output  is  segmfnt  poooe 

OOFiOOnOii' 
ocrjoont  i« 
00r,00nl i* 
OOFtOOftl  I* 
OOfiOOM  I* 

oor  tOOpi «« 
OCr^  OOnAi  5 
OOF ! OOpr I  1 
OOFiOOlf* 
OOr ,0D?5  t0 

oOF,ooiri« 

OOFtOO^C i  1 

ooeiOoi:  i5 
0OF|OOn*|3 
OUTPUTfOOE)  IS  005f  LONG 
2  00?i00o6«0 

00?|OOPA«0 
00?i00n6i0 
00?i0006l0 

2  oo?iOon«io 

OUTPUTSOLN  IS  SEGMfNT  0OO0F 
OOfiOOqO'O 
OOFtOOo""0 
OOr J00np|5 
00r J00o6|2 
00F:00n6l2 

3  OOFtOOpTiO 
OOF|0008|2 
00FjOOn9|2 

«  OOFiOOn«il 

oonooisn 

00r i  00 1 || 2 


«  oor  jOOi »  |2 

3  OOfiOOiCi* 

OOFiOO?AIB 
OUTPUlSOLNtOOFi  IS  003l  LONC 

2  00?i00o6'0 

00?iOOn*lO 
00?i 00O6I0 
OOpiOOfiAlO 
OO?700n6lO 
00?,00o»|1 
00? :OOft*lO 
00?t00nA lO 
00?iOOnA  10 

?  OO21OO06IO 

GOxpRYl  IS  SEGMENT  oofmc 

01 0 foopo  lO 
oioToonoiC 
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L*BTL  L 1 0  I 


LABtL  Ll.L?iL3iL<«L*«L».ST0Pti«l 


Ivi 


INTEGER  PRPcEOuRE  EUCl  I "»(  U.  Vl  I 

VALUE  U»V| 

INTEGER  U.VI 
OEGIN 

INTEGER  *• 

LABEL  L8»L9« 
Xl.ENTIE»(U/V)i 

ir  w»v  gtk  u  then  rEgin 

HI-W-1  | 
GO  TO  L<" 
FNOl 
IF  (**\ )»V  LEO  u  T"CN  »rGIN 

K  I  ■  ««  I  I 

fin    Tn    l'I 
Enpii 
FUCLIPIaHl 

ENn  EUCLlDl 


Ll  i 

«StOp  Ir 


I0TIMEl«TIME(3)l 
T!mi«timE(?ji 

PIV0TN(J"'«0I 

PIVOTNliMI«PIVOTNU"*l  I 
TIME  >  M»XTIME 
IF  (CUHTlHr*TlHE(2>*TlM1>MAXTlME  THEN  GO  T"  STOPI 
FOR  I •■»  STEP  t  UNTIL  H  nO 

IF  A(I.N)  LSS  0  THFN  RFf.IN 

fin  Tn  (.?, 
F»nl 
GO  TO  STOPI 
FOR  M»J  STEP  1  UNTIL  N«l  nfl 

IF  «{R.K]«0  TMFN  6"  Tn  I  «| 
GO  To  NUSOLliTIUNi 

NULLI«T«UEI 

FOR  Jl.K»J  STEP  1  UNTIL  N»1  HO 
IF  *Cs.J1<0  ThFN 
•  EGIN 

r  i. 01 

IF  *r  I.J'<At T.LI  THfN  Ll«J 

FLSF  Ir  »t  I.Jl«At 1,1]  THEN 
REGTN 

I  1*1*1  • 
GO  TO  L  3  I 
ENO| 
ENDl 
5  •  »0  • 

IF  A[S,L1«0  ThEN 
Rf  GIN 

Sl"5*l  I 
50  TO  LSI 


Lit 


LSI 


OlOtOOflOlO 

olo  ioono'0 

OlOjOOoOlO 

oioioonnio 

OlOjOOoOlO 

oioioonoio 

OlOjOOoOlO 

OlOsOOnO'0 

EUCLID  IS 

SEGMfNT  000 |  1 

01 1  ioonmo 

OH  tOOnoiO 

OH  i  0002  1 5 

on 'ioonti* 

Oil :00n*»0 

Oil t00o6>3 

0! 1  ,00oA  '  J 

Oil ,00001* 

on  ,0004 ,0 

01 1 ,0004 1 3 

01 1 iOOOA  13 

Oil lOOoRI 3 

UCLlOf Oil ) 

IS  OOOr  LONG 

OlOi 00001 0 

0  1  0  i  0  0  0  0  I  o 

OloiOOol  I* 

OtOiOOo3>2 

oiojooo*  10 

OlOrOOOA  10 

0l0|00Q4|2 

0i0i00n?i2 

Oioi00n8i3 

oin ToooQ 1 1 

010,00oni0 

010,00ori0 

oioioonr i3 

010lOOrtF«5 

010IOOOF»2 

0  1  0  t  0  0  1  0  1  0 

0 1 0  .'  0  0 1  «  i  3 

OlOj.00  1  6  «  1 

Oloiooir ii 

0  1  0 1 00 1 715 

o  i  o :  o  o  i  p ,  o 

Olo  t  00 1 o  i  « 

oi o  t ooi r  1 1 

OlOiOOir  i5 

0 1  0  i  0  0  ?  2 1  1 

Ol0;O0?*|3 

010:00?r  iO 

010:0090:2 

0 1  0  ;'  0  0  ■>  0  :  5 

oin  :Oo?e  is 

0l0i00?9>2 

010tOO?A|0 

010: 00?R I  * 

oioioo?r.n 

Olntoo?ni 3 
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.*  I 


17. 


Llm 


ENUI 

L*M8D*PEN0Mi«l I 

FOR  J I ■ J  STfP  I  UNTIL  L"l • L*1  STEP  1  UnTU  N-l  nO 
IF  »[»fjl<0  THEN 
•  EGlN 

FOR  I  l  -0  STEP  I  iiNTjl  S-l  00 

If  »fI.Jl  Nf8  o  THFN  GO  TO  L7I 

Tl»rUCLIf)(»[t»J)f»(SfLl)l 

ir  r  T«»rs»Li"*f  s»  Ji )  ano  (t>d  then 

REGTN 

T  LSI 
T  t  •  T  *  1  I 

ir  T.»t!,LiP»ti»ji 
then  r,n  to  L* 
ELSE  IE  T*At I • L  3  > A  C I. Jl 
THFN  Ti.T«l | 
ENni 
IF  "*tRi Jl*L*MR04DFN0M  >  LAHRO*NUM«T  THF.N 
RE6TN 

L»MRO»NUM|.-»f R. Jj  I 
1*M«0»0EN0M|»T| 
ENOI 
END 

El  SE 
MULLI.F/kLSEI 
C'>rUCLIO(ltR«N)»LAHB0*nENOM.L4MRrtNUM) i 

s«"-(c*i )» 

T • "-At  R  »N  J  t 

»(FL*6I»TRUE| 

IE  MULL  THEN  GO  TO  HOI 
FOR  Jl»l  STE^  1  UNTIL  N»l  00 
IE  *CR.J3>« 

THEN  .EgIN 

Cl«EtlCLID(»tP. jl*L»HRO»nFNOM.LAMR0ANuM)| 
IF  S.At««J]<T»C 
THEN  REG'N 

TI>A(RtJll 
Si. 'I 

NFL8GI«F«LsEl 
FNfM 
ENOI 

jF  ^.LAMBD*KUM  <  TtLAHROAnEN'1*! 
THEN  rEGIN 

LAMBD»NUH|.IF  NFLAr.  ThEw  I00»y-1  F"L  jE  T| 

LAMBOAOENOHI-IF  NFLAG  THEN  10o*S  ELSE  SI 

I  NO  I 
EOR  Ji.l  STEP  1  UNTIL  L  - 1 » L* 1  STE"  I  UNTIL  U    Op 

BEGIN 

r i ■  t r  LAMRnApfNOM  nEo  o 

THEN  EUCL!D<»rPiJl*LAMBnAOENnM.LAMBOANUM> 
ELSE  IE  »tR.  il  «  0 

ThFN  -1 

El«F  Oi 
IF  C  NEO  0  TH?N 

FOR  Ii.O  STfP  J  UNTIL  N  no 


OlOiOO-JF 

oio, oo?r 
oio  ,oo*o 

OlOtOOlO 
010,  OO^n 

010  ,001F 
OlOiOOlF 
OlOiOOat 
010,00«* 
OlO  ,00a« 
01(1  i  00a  F 
OlO i OOoF 
010:00*0 
010,00*? 
0l0:00^S 
OlO|00*S 
010,00*9 
010 ,00*R 
01 0  jOO«n 
Olo  toosr 
Oio  :00*r 
OlO, 00*0 
010,00c! 
010,00*1 
OlO  :00*1 
0l0l00*1 
OlOlOOnS 
OlO I  00*9 
01 Oi 00*  i 
OlO, 00*r 
OlO  !  00*p 
OlO I 00*n 
010iOO*F 
0l0,O0*r 
OlO ,00  TO 
OlOlOOTO 
010,0076 
010:0077 
010,0079 
010^0071 
OlO -007P 
0  1  0  :  0  0  7  R 
010:00  7" 
010,0070 
010  .00  7 r 
OlO J  00*0 

oioToopi 

01 01 00«* 
0  1  0  :  0  0  M. 
OlOjOOol 

oio^oor  < 
0101OO93 
0107009* 

0 1 0 ,"  0  0  9  9 

0  1  0  .'  0  0  O  A 

010 :00oa 

oioTooof 
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M  1 .  j]  i.«t  r .  J!»c»»r  !.U  > 

EN0| 
GO  To  111 
STORl 

T!«i.(T!«<E(2>«T!M)  OIV  «0| 
inTtHf l.( TlKt(3)-I0TlHt)  olv  »o, 

end  soxnuYi i 


llltltlt 


STE 


no.;OlmTi 
*E«T, 

MtlltM 

STn.Pi 
CNn.  i 


tiiiiiiiiuiittuniniii'iititKiiiiiiitituinitiitiiiiiiiiKi 

llltt1tlllIllIlllt!II!lll<llt)!)lll>lllllMIIIIItllIltlHItlII|t 

CUMTIKf 1.0  P 

CRE»TE.OJl 
•1 { ! TER. 1 . MU«ITC«  >  RtBIN 

TI LLMATHI K | 

OUTRljTi 

GOMORTj | 

OUTPuTSOLNl 

■»I TC< L 1"E  > I 

HRITEUJUEt^PRoctSSOR    Mmc    in    SECnuns   «•!•». Timii 

*R!TE(LINE,<"!/n    TIT     I"    SrrnN«^».IIO».!OTI><E)l 

N«ITE(LiNE.<«NUMBE»    Or    MVOTt    "  .  I  »>  .P  I  VOTNIJM  )  I 

IF  CUMTIHEl«CUMTIME4TI»*»0  r,ro  M»»TTME  THEM  GO  TO  STOPl 

GO  TO  hEXTI 
INI 

KR|TE<lIKE.<"NO    SOL|iriQ«">ll 

FNOI 

'timninttiimiiaiiimiiiiiinuiitiiHimiimtiiinim 


number  Of  errohs  orTri-TE''  •  aooo. 

number  or  segments  ■  ou.     tot»l  segment  $ur  •  cor?*  »n»os.     cn»t  estim»te  •  oo?J3o  nOros. 

PRflGB*"  SI7E  »  000S09  C»»OS-  0025'S  S»NT»fTIC  ITE'S.  OnOO««  PIS*  SEGMENTS. 

pROgR'm  r  I  i.  E  nImEi  0ptIm"l/»RE' . 

COMPILATION  THE  •  0000|«  SEfONDs  EL»PSCO.  00000*  SECONDS  PROffSMNG. 


010 |00«ni 3 
OlOiOOlSl 1 
)  010(OOi»|3 

01 r , 00i9i 0 
0!0i00i9t0 
010. 00. „, « 
O1o':00aEi2 
G0"PRT1 (OlOi  IS  OORl  LONG 

?  oo?ioo«i»m 

00? i00n6 10 
00?fOOftt  10 
00?j00n*l« 

00?  :Onnr  i? 
!  00?,oon»iO 

00? lOOndiR 
00?iO0n9l2 
00?|00P»  10 
00?|000«'« 
00?  i00r,r  l  1 

oojtootnil 

00?,00?«i  1 
00? ,00l«l  1 

oo? lOotr i2 
oo?  iooir is 

00? |00*C i' 

00?,0P<?|2 

00?,G0«?|2 

2  00?iOOiai< 

00?  tOOat ■* 

00?  |00«« I* 

BLOCK(00?)  IS  OOAn  LONG 

0»T»  IS  00?^  LONG 

ST»fK  f sTIm«TE.O0P««S 
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APPENDIX  B 
HEURISTIC  AREA  LOCATION  CODE 


ioU 


ijRRnur.HS    i-ftrnn    »i  r.ni    cnM"!tro    lfvEI.    1 1  ,n?.O0,0l  s. 


fRIOay,  O6/30/??.  (1"Ji3a  am. 


*  r  r  a  / 


EnRiSTlr/ntf 


E  f  T 


pEt.'n 

•HILE  TRUE  pu 
••  E  r, »  n 

t  TUNING  o|R»mETFBS  * 

nEFTNF 


MAYMi|MP»RrFLS" 


t, 

t    THE    maximum    miMPfR    nr    PARCELS    A| L n«En-uSFD    Tn 

t   LiEminsIdn   arrays 

mAXMumCNSTRMS" 

n 
f . 

»     THE     MAXIMUM     NUMp'R     nr     CONSTRAINTS 
uAx«i|MMqnRS» 
u 
f  . 
f     'HE    MAylMUM    NijMpFR    nr    AO  /ATFNT    ►>ElGHRtl«S    THAT     A 
I     PARTEL    CAN    MAVE 
MA»Mi|MP|ir)[)tES"MAXNIIMpARCELS     *IV    ?    $1 

1    THE    MAXJMIJM    wiiMPrR    nr    pimpLFS    EvPECTEO 
!«ll«ttt'itUnMUI>tllMttMU(<IItMltU«l|lltllU<lt>IIUI<Mlltt>t> 

iEF'nF  LlNE»iTATinN  fl 

rlLr    ST  AT  inK(Kl*Jn, REMOTE  ,MAx»ECST7f  ,12.MV"SF.I0)| 

rlLr    nIs«<MNO«"IS><.MAXPFCSI>F«1/i. BLOCKS!  ZE"«20>I 

rnR"»T    F](U>X},iJly,,i?tX9in9,1?)l 

lNTrGER    ARRA»    T-Pl  r  1  |M4xniimp»rcfi  <;! . 

TMp2r1 |MAxmump«RCELS1. 

TMP^ri|M»XNIH1C,'STRNTS»tlMAl'NUMP4RCFLSl' 

TMPItltM'XNllMpARCEl.S.llMAXNllMNt'nBSI. 

TMPSri|MAXNIjMrvSTRNTSl| 
nEr  'sir 
■  Tf  BAoCELKUf-  !• 

T»P1fPtRCFLNIIM'      I, 

»  tn   array  ihat  holos  the  iNTrr.ER  weights  inoEyfu 

»     *Y     TME     I'TFPMAL     DAHrr|      MIlHOFRS 
pUO^LFNuMf  PJKCEUNIJMlr 

TMP?f  PA  *CELMIM1     f. 

j    tn   sHpAv   holding   cuporiuT   PunniE    in   num^Eu.    ir   value 

*     IS     7FRn    THEN    PAR'FL    nnFS    NHT    RClPNG    TO    A     PUOni  f 
rNSTRNTrNI,M,RARfFLNIIMJ« 

TMPlf  vUM,p*RrELklliM]     l, 
»     H    *gpAY    HO|MNG     THF    rnNsTR«tM    rnEFFlriFNTs 


000:00n0t 0 

RLHCt  IS  st-GMFNT  000n2 

1 

oo?looooiO 

00? :0000<0 

? 

oo?.oonpiO 

00?: OOnni 0 

00?;00o0«0 

00?:00n0»0 

RLPC"  IS  SEGMFNT  00003 

ooi; oonoi  o 

003:00001° 

OOJiOOnOtO 

003 ,00o0t 0 

003,0000(0 

003 ; OOnO i 0 

003:000010 

003:0000' 0 

003 |00o0|0 

003t00o0|0 

003:000010 

003:00o0>0 

003.00oO|0 

O03't00o0|0 

003;0000«0 

OO3tO0oOt0 

003  lOOoo  lO 

003iOOoO»0 

003,00o0|0 

DATA  IS  000*  LONG 

003i00o0l0 

niTa  IS  OOOa  LOnc, 

003t00oOl0 

003tOOoOiO 

003:'00n?|0 

O03tO0n«|O 

003:000711 

003;00nii2 

003-OOnr |2 

003.00ftf!2 

O03:00nri2 

003,OOOC;2 

003:00nr :2 

003i00or |2 

003,00or:2 

O03,000r|2 

00i:00nc  <2 

0O3t0O0C«2 

003;00op|2 

003,00nC(2 

105 


t 
t 

RHSrNUMl» 

TMP 
X 
I 
I 

NBO'fNUM.PAR 

I 

X 

INTEGER  COUN 
I 

INTEGER  NUMO 
I 
I 
I 
I 
I 

INTEGER  ARRA 
t 
X 
I 

INTEGER  "ATE 
t 

INTEGER  NATE 
X 

INTEGER  ARRA 
t 
I 
I 

INTEGER  PDLI 
I 

INTEGER  NEXT 
X 
t 

INTEGER  ARRA 
» 
I 

INTTQER  SOLN 
X 

INTEGER  ARRA 

f>Er!NE    CNSTR 
t 
I 
INTE'.ER    array 

nEr»NE  cSUMt 

* 
* 
% 

!  ARrL  STnpi 
uOOLEAN  SOLN 
t 
* 
nEFTNE  STFP1 
»  *  (  X »  %  X  %  »  *  t  %  t  X  *  *  I 


FOR  EACH  PARCFL.  THFPF  IS  *  CtlEF,  FOR  EACH 

CONSTRAINT  (FIRS'  INDFX)  OVER  EACH  PARCEL  (SECOND  INDX) 


5[NUMl  f 
A  AdrAy 
RIGHT  HA 
CESS  OR 
CELNUMJ. 

AN  ARRAY 

NEIGHBOR 

Tl 

CURRENT 

FPARCELS 

RESPECTI 

NUMBER  0 

NUMBER  0 

PROCESSO 

INcLUDE 

Y  NATERL 
AN  ARRAY 
SMALLEST 
DEPTHS  0 
RLEVELSP 
POSITION 
RTABLEl 
CURRENT 

Y  POLIDN 
CURRENT 
THE  INDE 
HAS  NO  S 
DNUMPOSI 
POSITION 
PUOOLENU 
THE  LAST 
OR  C0M8I 

Y  SOLNSt 
THE  PUDD 
MEET  ALL 
SPOSI 
POSITION 

Y  Tmp7[0 
NTNAMEC  I 
THE  T2  r 
THE  CONS 

TMP6I » l 
POLIDNUM 


INDExED  4Y  COMSTR* 
ND  SIDE  v'ALUFs  THA 
GREATER  'HAN 

TMP4fPARrELNIIM,NUM 

HOLDING  THE  INTEr 

S  OF  PARCEL  PARCEL 


I nt  NumoEr  That  holds  ThE 

T  FACH  CONSTRAINT  MUST  BE 


1  tl 

NAL  PARCEL  NUmrErs  OE  ThE 

Hum  (SECOND  INDEX) 


LEVELS  ARRAY 
I  HE  J 


RFDi 


INDEX  TO  THE  WATER 

.NUMOFCN^TRNTs.PRT 

VELY 

F  PARCEL*  TOTAL. 

F  CONSTRMNTS  ENTE 

R  TIME  USED  Tp  LDC 

nATA  INITIALI?ATI0 

EVELStll"AXNUHPARC 

HOLDING  UNIOUE  WE 

TO  LARGEST.  THF  V 
F  THF  HATER  LFVELS 
OSI 

OF  THE  L*ST  FNTRY  IN  ThE  NATERLFVELS  ARRAY 

WATER  OE»TH 

UM[1 IMAXMUMPUODLESI r 

ID  NUMPE»  Of  FACH  PUDDLE. 

X  INTn  THE  ARRAY 

I6NIFICAYCE 


ATE  A  SOLUTION  (DOES  NOT 

N  TIMF) 

ELSII 

I6HTS  ORDERED  FROM 

ALUES  ARE  THE  SUCCESSIVE 


IN  THE  POLIDNUM  ARRAY 
ASSUMED  TO  A  NEW 


OF  THE  LAST  ENTRY 
Ml 

PUDDLE  TO  NUMBER 
NED  PUODLE 
llMAXNUM'UODLESI I 

le  m  nu«rErs  of  the  puddles  that 

THE  CONSTRAINTS 


OF  THE  LAST  FNTRY 
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005t00?s t« 

005,00?A  |« 

005,00?ai« 

00?,00?AI* 

005,00?A|* 

005(00?AI« 

005,00?AI« 

005,00?A»» 

A 

00«.tOO?AI« 

r.FTNBHRS  IS 

SEGMENT  OOOoD 

OOmOOOOlO 

OOn  j  OOOA  t  5 

00nt0004  |5 

OOn  t00O6l3 

00nt0006i 3 

00ot00l6|5 

000:00l''i3 

OOpjOOi Tt  3 

OOP,OOlClO 

00o,00?5|3 

00ni00?5|3 

OOotOO?AlO 

5 

OOn  |00?E'3 

0OPiOO?Fll 

00n,00l?|0 

00P:003J'«5 

00p  tOOin  i» 

OOP  ,00a  Ji  3 

00Pl004fl|2 

5 

00r .OOos i 2 

GFTNpnRS(OOD) 

IS  0050  LONG 

4 

005:00?At« 

005T00?Al« 

oo5,oo?ni2 

005,00?ri0 

005 iO0?r i« 

TlALlZEnATA( 005> 

IS  003n  LONG 

3 

00Vt00l7i3 

003:001713 

00lt00l Tt 3 

00?i00l7«  3 

3 

003:00)713 

OliTPllTwTS  IS 

SEGMENT  OOOOF 

OOf : 0000  !  0 

OOr :0001 t« 

00F:00n2i2 

OOr ,OOo7i5 

OOF .'00n7i5 

4 

OOFiOOnCl? 
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STFPHJ.l.N)     WTrTFMP|«TFMP*1])l 
FNn  i 
J  uEm<"mREk    T  I  Kt     *"EN     START     LOfl'INr,    FOR    Sni.nT  I  n  m  S 

PRT''1FI  «TJH(  (  ?)  I 

fno   nil  I  put  w  t  s » 

oRtTCnURF    NF»TK»TFRT»RLF| 

*  SUPPLIES  THF  *'FXT  *ATrPTABLF 
•EG'N 
»      fHF'K  FnR  Nn  SnLUTION 

|F     "01INTl«CPUNT»l     r,T9    «A  T  FRL  E  VEL  «POS 
TMFN    BfulN 

uPMFf  I  TNF,<//t"Nn    SOli'Tlnw">)| 
PRlIMf  l»(TIHF{?)-P'TI»»F)     *MV    ftOI 
*BJTE(LlNE.<//."PR''lCLSSfiP    TI»«E|"»T5»"    SECS,">.PPT)ME) I 

gh  Tn  sfnpi 

ENIH 
«»T('RT»RLEI«»4T''RLEVFLSfrnilMTl| 
FNO    NEKTwATEKTArlEi 

oRn*Eni)RE   l.  CATFPunntFsi 

>  uetfbnins  which  »»Hrri s  FtnnnFO  and 

»     AGGREGATES    THFM     iNTn    PliniUEs     fl.F,     AREAS) 
"EGT»j 
I AS^l    NFXTPAKFEL I 


OUTPUT«TS(OOF) 
3 


INTEGER    I.  J.K.PARCFL.TEMP1  I 
TNTFGER     ARRAT      TFHpr  1  I  M»  X  NIlMN  "QRS  1  , 


I nCATFPUDnLES     IS 


pRfl^EOURF     NF"PU«nLFi 

*    C»E»TES    A    NTV.    PUT1LF    MAnE    Or    A    SINGIF    PARCEL 
'EG'N 
oUO"»LENll«rPANCELl!"MFXTPIiOniENUM|r,NE,<TP(innLFNIIM*1l 

oOltonum(pi>i  lnN>iMpOs,«pruinN  iMpp«4.i),«NExTpur)0LENu><, 

rNO    NFWPIIIinLS  J 


oRo'FnURE    CtiHti  [MFpUnnLES| 

*  ccimrinEs   the   uEwtr   Fi.nOr»En  parcFi    with   all    TmF 

>     LXJSTINr,    PIJf'OLFS    THAT     NFTGHqflR    THE    PARCEL     ANn 

*  OIVr<;    THF     E^PsNflE"    PiinnLF     A    NEW     ID    NUMb^R 

-EG'm 

anrii.FAN    FlRyTl 

tEmooRArIlY    ASSIGN    NF^LT    FLmnEn    PARCEL     TCI    fiNE    Or    ITS 

ME  I  ".HRORlNG    ►'MOHLES 

oUmLENuM[PA«CELl'"''FMPM  1  I 

ASST'.N    ALI      ThF     PARCEl  S    PFtnNMNr.    TO    THE     PIIOnLES    THAT 

aur    jt\    of    Cl'MHlNEn    TmE    NFh    1  ">    NiimhFR 

k'EXToijunl  >"Ni.it,MFXTPnnni  fniim*  l » 

<;  T  F  «  1  f  I.  I  ,nI'«(UfpaRCFI  S) 
STFPKK  •  1  ..I) 

If     PlinnLFMiHf  T  )    EOL     Tfmp(k1 

THEN     piinnLFNH"!  IllsNrXTPlinnLrMJHI 
uPD»TF    pni.Il'NUM    ARRA*    OY    GlV'NG    THE     to    MIMRFPS    OF     Tht 

rnn3jNFn    PutwlFs    thf    nF«   Pimv.E    »>i|MRFR 
«TF°) (I.I.J) 


CnMRiNEPunnLEs  IS 


OOr  tOOOF lO 
OOfiOO?! •• 
OOF|00??l 1 
OOFtOO??l 1 
OOF :O0?3»5 
IS  00?*  LONG 

001,001 r  1 3 

003,001  n3 
oostooi ri3 

O03i00)7l3 
003i00iri3 
003,00)7«3 
003|0019|2 
001,00) Al2 
OOl,00)Fl5 
003,00??, 3 
O03,0030l« 
003|003?|2 
003i0032l2 
003:00l*i0 
003l00l7l« 
003i0037l« 
003,O0j7i« 
003|0017»« 
003,00371* 
003,0037i« 
SEGMENT  000)0 
OloiOOoO'O 
OlOjOOflOsO 
OtOiOOol  ' 
OlOiOOol :' 
OlOiOOnl : 

oiotooni  »' 

OtOlOOOl I 
010.00oS,0 

01o|ooo7|S 

OlntOOpSiO 

oi n ,oons  |0 
oio,oon*tO 

OlOjOOnAiO 
OlOiOOOMO 
OlOjOOoAsO 
OloiOOOAtO 
SEGmfnt  OOOjl 

OH  ."OOnoiO 

OJ 1 iOOOOlO 
01 1 tOOflOiO 
01) rOOo?lO 
01 )  :0Op?iO 
01  I :00«?lO 
Oil (00o3l2 
01 1 . 00n3i 2 
on  t00nri5 

Oil  lOOOF  10 
01)  lOOt  1  '5 
01)  ;00l  1  '5 
01 1  [001  1  i5 
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sTfpi  (  k  •  1  ,"ni  inniiMPns) 

IF    POLI^NiiMfK)    F«L    TEMPrIl 

ThFN    pOLIPNUM'Kli«NE)(TpUnPLENUHI 
rLE«N    UP    POllilNiiM    4PP4Y    HY    F  tIHTM»TIMfi    PiiPLTCATE    PI'PPLENiim  '  s 

K   I  »  '    I 

r  imri'TPijEi 

sTE°n  I.  l.PDLIOMUMPn<;) 

PET,  IN 

IF  PDLIUNU*(I]  NFO  NFXT»UOP|  FNHM  np  FIpsT 
THEN  BEGIN 

if  PnLinMiHtn  Eni   nextphddlfnuh  then  firsti.falsf  i 
pnLlnNiiM(Kii«»nLinNiiMrili 

X l«K»,  I 
ENO 
El  SE  J i ■ J+  1  I 

END| 
pDL'DNUMPfisi.PDLIONMMPns-JI 
FNrt  cnHnlNEPUDnLESI 


COMPTNEPUUPl  ESf 0! 1 ) 


STE°1(P»RCEL.J .  NUHOFPARCELS) 
P  E  r,  I N 
i  FIND  NF«  PARCELS  UNDER  "ATrP  « 

IF  «T[PARCELl  *F1  WATER'AfJLF  THEN  r,D  To  NEXTPAPfELI 
»  CHFCK  IF  NFIGHrORS  bELIKT,  TD  punOLFS 

»  save  the  puddle   id  nps.   if  the*  on 

Jlal  I 

STFP1(I»l#M*XNllMNRnRS) 

IF     TE"P1 1-NflDRt I .P'RCEl 1     GTR    0    THFN 
IF     TEMPt  .11  l»PUOni.EN|iMrTFMpt  1     GTP    0 
THEN    oFGIN 

*  CMFCK     IF    ID    NiiHRE"    MS    SAVED    FJEFPRE  5 
KJ  m\  I 

WHILE    TE^tM    NF<J    TFMPIJl    00    K  I  aH  •  1  I 
IF    k    FqL    J    TnrN    j i > j«  1  j 
FNOl 
IF    Ji-J-1     FQL    0  5 

i  NfjGHonRS  nnNT  belong  m  puOPlfs.sft  up  ht*   pupplf 

THEN    nE-JPIIPOLF 

♦  CpMrtlKif    ALL    PUPPLF*    TH»T    nFIGhBOR    PARCEL 

Else   rnMRiMEPiionLF,;i 

mF  <TP4KCrL I 

FNn  I 
rMn     LHCATf  PlIJULFSI  « 

LnrATEPUUPLES( PlOl 

°nn*rnupF  fHtc«CNSTwwTS: 

t    >n«    F»f.H    >>hppi_F    'HE    PpprFDUBF    ChFCk    IF    M  L    ThF 
»    ii)N?TO«InTS    ABE     "FET.       TF    PnF    Op    MORF    PUPDLEs    SATISFY 
?    ,,ONSTRA  InTS    TmE>i    A    SfUljTTIlM    HAS    pECN    FDiiND    Ann    procfss- 
»      ING     STPP^ 

<F.r,i  v 

|AH"i.NFXTPlMPLe':  5 

chECc nstpnTs    IS 

lNTr,FR    "ARC  I  L.T.J, k , 

FUR    UC»    rflrSTWAIMT     Finn     THF     SUH    pr    THF     rPNSTOAIwT     rpEFICTfNTS 
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01  1 
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« 
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on 
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2 

on 

00?6 

5 

on 
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2 

on 

00?9 

2 

on 

00?r 

0 

on 

00?F 

3 

on 

00»F 

5 

on 

00?F 

5 

on 

003t 

4 

on 

003? 

1 

on 

0033 

4 
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OOPA 

0 
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rOOoR 

0 
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OOPS 

0 
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OOor 

3 
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3 
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OOnr 

3 
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OOpF 

3 
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OOOF 
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OOPF 
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1 
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OlP 
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|3 
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2 
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10016 

i2 
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r  00  1  7 

0 
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lOOlp 

3 
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[OOlF 

1 

010 

I  00?0 

i3 

010 

|00?1 

,5 

010 

:00?1 

b 

010 

■0093 

il 

Oln 

:00?3 

tl 

Oln 

|00?« 

|3 

Oln 

:00?4 

i3 

010 

•00?? 

■  0 

IS 

i0?n 

.Owe, 

003 

:003f 

|4 

003 

:omr 

14 

003 

r003T 

■  4 

003 

:003T 

|4 
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•  00i7 

|4 

003 

:0037 

|4 

003 

:00  3  7 

:  4 
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;0037 

|4 
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FNT  P 

D0l2 
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lO 
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:OPnO 
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»         of    r»cH  parcel   tn  eath  pnnnir 
«Tr°i  r  k.  i  ,PnLin*jiiMPnsi 

STrpl(  I.l.MljMnrrNSTPNTSl 

/•Suit  k  ,  1 )  i  »oi 

<;TE°1  (PARCEL  •LMllMflFpARCFLS) 

IF    JiaP  jODLtNtJM[  PAPCFL1    GTR    n 
ThEN    "ECI«J 
<  l"1  I 

while   J  nfo  p^LinmiMrKi   nn  ki»k*ii 

STEP  If  I, 1 ,  NUmVCNSTRNTSi 

CSUMrK.Ml«CsilMrK.IJ*CNsT(»NT[  I  "PARCEL  II 
FNPI 
t  rHF*<     IF    ANY    PlinOLfS    MEET    »L.l    CONSTRAINTS 

<;0L,lSPnsi.0l 

<;TE»ifK,i.pptinwuMPns> 
*f  c.  In 
STEPK|»l»>IIIMfireNSTRMTSl 

IF     CSUMIK.T1     GTR    R-SIM 

Then   an   to   NF'TPimnLn 

SOLNibTSUEI 

soiNst sulnsposi«SOlnspo«*ii i »pnL IOnumik i i 

POL  IDNIiMfK]  (  »  0  J 
•JExTpuOOLEi 

FNOI 
FNO  CHFCfCNSIRNTSI 

KiitiiKiiinuniintttnuttts'UtiiituiiutMUiutittuiiMiiiH 

pRO^EOURF  OuTPUTSOLNSI 

»  PRINTS  OUT  SOLUTIONS  ANO  OTwfP 

*  PUOnifS  THAT  FAILED  Tn  "EET  ALL  THE  CONSTRAINTS 

oEGf  M 

lNTrr,E«  PaRCELi 

1RR«»  A  rO  11 1  ) I 
Dni«lE«  PTRAi 
!NTr(iEp  NiSLi"'!*  J>K| 


PRO'EOIIRF     WRlTAf  VALII,  IfPTB)  I 

%   manhlES  variable  format  nr  output  line 

V*L    IE      ■'Al  U.  I  ! 
REAL    VALU| 
lNTrf,ER     II 
of)  I  M  TCP    PTRi 
>  F  C,  '  N 
S"ITCM    FHRMAI     FT  I .( A3  ), ( I  3 ) I 

vAL  'E     ARPAY     ETSlZEf  3.3) I 


.MI T£(PTR.F  T(  I  )  . VALU)  I 
eTBi.PTR*f  TSI2E r  I )  | 
rM'i     «RI TA  ; 


iifw-MRER    PROCESSOR    T  t  ME     FOR     r|NOT»'G    THF     SOLUTION 
pmt 'MF '«< TI*t <2 J-PRTtmf )     njv    AOl 
*»IH     ALL     IK     PmniFS 


01? 

01? 

OJ? 

01? 

01  7 

01? 

012 

a 

01? 
01? 
012 
012 
01? 

« 

01? 
01? 
01? 
012 

a 

012 
012 
01? 
012 
012 
012 
012 
012 

* 

012 

CmfCKCNSTrntS(012) 

IS 

X     3 

003 
003 
003 
003 
003 

3 

003 

outoutsolns 

IS 

SEGm 

013 

Oil 

013 

013 

013 

013 

013 

013 

013 

013 

013 
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013 

-RITA 

IS 

SEGm 
014 

OATa 

IS 

RATA 

IS 
014 
014 
014 

*RTTA(01«1 

IS 

It 

013 
013 
013 
013 

lOOnOlO 
lOOon  iO 

tOOpOlO 
i 00n«  i  3 
:00firi2 
tOOpc i2 

i 001 21 3 
iOOi JI2 
i  00 i a  i  0 

,ooirt5 

tOOlTiS 
•"00?2l« 
iOO?3ll 
tOO?3i 1 
|00?3<> 
|00?S|» 
l00?8l2 
,00?B|2 
i00?ril 
|0030|« 
tOOM |2 
i0034 iS 
t0036|2 
t0036t2 
1OO361S 
003a  LOnG 

tooiri* 

|003TI« 

iooiri« 

t00l7l« 
100371* 
l003fl« 
FNT  00013 

toono 10 

|00o2|0 

[OOn?iO 

|00n2l0 

1000210 

i00n2i0 

:00o2l0 

t00n?t0 

iO0o2«0 

iOOn?io 

:OOo2lO 

!0002I0 

ENT  00014 

:00o5l3 

0051  LONC 

POO?  LOnc, 

iO0oSi3 

|0013«3 
i00iSt5 

oojr  long 

i00o2«0 
(00o2|0 

i00n2»0 
i00n4  I* 
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NI»«8RT(NUMPF PARCELS) I 

uRI  TE(STATlriw<//»"»LL    THE    A»EAs    ARFI".//>)I 
STEP1(K.1.N) 
BEGIN 

PT«A|«P0INTER(A )l  « 

STEPKJ.l.N) 
"EGIN 

PARCEL l»2«*(**l )*J|  5 

!r    PllonLEN|)M(»ARCFL  ]    GTR    p 
THEN    rEGTN 

WRJTACwTf PARCElI.I.PTRA) I  6 

FN" 
ELSE    «RITA(-    .     ».O.PTRA)l  f> 

ENm 

WRlTEfSTATION.l?.Af*J )l  5 

ENQl 
WRITE  ALL  THE  SOLUTIflNS  < 

wRITE(STAT10M»<//»"AREAS  THAT  MEET  CONSTRAINTS  AREI»,//>)| 
STE»HI.1.S0LNSP0S) 
RESIN 

WRITE(STATION)l  4 

SUMI.OI 
STEPl(K.l.N) 
BE&IN 

PTHAIaPOlNTER(A)!  5 

STEPl( J.I.N) 
BEGIN 

PARCEL<"2"*<K'1 >*J»  * 

IF  PuODLENIlMf  PARCFl  j  E8L  SOLNStI] 
THEN  REGTN 

NRITA(MTrPARCELlt l.PTRA)!  7 

SUMi«SlJH  +  wTtP»RCELu 
END 
ELSE  KRITAC"  1  ".O.PTRA)l  7 

ENOl 
WRlTEf  STATION.  12.  *  f  *  1  -> »  6 

ENOl 
WRlTE(sTATTON.</.«sUM  OF  WETp.HTs  ■  «•  1 1  0>»  SUM)  I  5 

ENni 

write  out  puddles  that  arf  n*t  solution  « 

»nd  the  constraints  that  the"  eatl 

writef station. <//."typf  yes  'f  you  want  to  knoh  which  "./. 

-CONSTRAINTS  Are  vIOL*TEn  rnP  ArEAs  NOT  IN  THE  SOLUTION,-./, 
"TYPf  MO  0THfRi«ISf'"«/'"I".«"O*">'I 
rEA^C STATION. 12. AC*  1)1 
jE  BOINTEP( A)»"yES  " 

TMEN  rEgIN 
wRITEf STATION. <//, "AREAS  THAT  vimATF  CONSTRAINTS  *RF I ". //> ) I  « 

STEB1(I.1.PPLI0NUMP0S) 

IF  P|)LlUNUurn  GTR  P 
THLN  qEGlN 

SIIMI.or  5 

KPlTFtSTATinN''  I 
STEPl (K.1 .Mi 
"EG  IN 

PTRA laPn'NTFRf A  > I  6 

STEP1 ( J.' .N) 


Ol3iOOoal « 
Ol3i00p«l2 
Ol3|OOPRt« 
01 3:00pri* 
01 3  I  00 1 C  l  1 
01 3 lOOl  1  '* 
013,001  1  l« 
01 3  t 001 6  1 1 
01 3 i 001 8  I  3 
013,0019|3 
013:001 Al2 
0l3«00lC'« 
013I001C* 
013:001FI5 

01  3  :00?0t2 
013|00?5«1 
0l3|00?5l« 
013|00?5«* 
0l3:00?Bll 
0l3t00?RI  1 
013l00?FI« 
Ol3l003««l 
01 3 i003«l5 
013i003«l5 

013tO03««2 
0l3«003Al5 
0l3i003At5 
0l3,0C-lFl2 
Ol3tOO«l •« 
Ol3lOO«3l2 
013tOO««l2 
0l3:00«6i« 
013jOO«R«5 
0l3»00»8t5 
Ol3tOO«Bi5 
01 3t00ACl2 
01 3 t 00^1 1  1 
013,00ms* 
Ol3:00?FI« 
013:00*0:1 
013:00*01 1 
013:00*0:1 
013,00*2,0 

01  ■<:  00*2:0 
013:00*?«« 
013:00*R,0 
013:00*ctO 

oi  3  ioo^ri i 

013:0074:* 
01 3: 007  4  |« 
0l3,'00r4:l 
013;007flt0 
013:00tr:« 
Ol3!00»o:l 
013:00*0:1 
013:00*4:4 
0l3-0O«6:  1 
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nff,  IN 
ri»-'lli'('(«-| )«J| 

ir    oiinm  rn'iMrfjprEi  i   i«i    pplIdnhmIji 
t«f.     >fc.t» 

J»ITt(»Tr P1BCFL 1.1  .PTQ1  )  I 

ni«nn»»»TtPiiirri  )i 
fmo 
ti  tr   <b|  t »( "    .   ",n.P79»  >  i 

FN"  l 
•  ipri  STi'Inm,  t?.M  •  l !  i 

•  j!  tf  (  sT»T!nN-«//.»r nNS'oSliTs   Th»T   irE    vlniiTCD" 

-   »ori">)i 
*TrPi»«.|.>'UMire«i<TP«TS> 

IF    rSuHf'.id    pTP,    «"Sf«? 

tmF>    «PiTF<SMT!nN, <!?»»>. 

STFPWJ.1,1?)     TMPfr  12*(K«1  )«J-1  I 


II 


■■»!' 


Ktttttt' 

TM' 


rxni 
,-PPPrrsSOP    TJT    .    -,io. 


<ce«.">.ppu»E  >  i 


Ult 
MT«T 


ituiinniiiiiivi'iiiitiinniimtini'iilliitimiim 

nilli 


f    NOT    set"    "'I 
J£MK 

Nf'T 

LTr  J 
phfc 


I1TI>  TiHLf  I 

rEPuOOirsi 
■  Ci-iTPMTSI 


'IjTSnLKSI 

iiiiiiiii>iiimiiiiiiiiiiniiiiimiitii)iimitimimiiiiti< 


NUMPfH    pi    fapnjs    iiFTr'Tfi    ,     .-ion. 

numbE*   pi    sf'.-IN'S    t   ii? I.      T>iTi.i_    sli.-rxT    <t?'    =   nMli«    .ions,      cut    f<.T1»»te    „   dpso'p    »p.pps. 

HPGB»M    S 1  ?  f     •    (-•"(, sM>    CU'i^i     PP?>;>?    (iMIrlK     T  Tr  •"; .     TO"*?    PT*«    tfr,»FMTs. 

*K0r,P*M    flit     M>ri      »pF  »/ jfrPUTI,  /Hh  if  f  T. 

ttlMPIlM  ll'K    1  I  it     s    i)r.P(L*    <rnui.l    tL'i'Sfl,     iv.ppp,     SfPPMPS    o»n(-FS<  1  *C. . 


01  1 

00«»: 1 

r 

013 

OOPll* 

013 

OO.mO 

013 

00«r  i* 

f> 

Oil 

O0«FI* 

111  3 

00»?|0 

"13 

OOOlt | 

* 

013 

PO»«i  1 

01 3i00«T  |5 

7 

01  3:00am? 
01 3  j  OOpni 1 

A 

nl3,00«pt«      . 
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APPENDIX  C 
NARIS  SYNTAX 
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1.   Fundamental  Definitions 
<empty>: :=  a  null  string 
<letter>::=  A | B | ... | Z 
<digit>::=  o|l|...|9 

<special  character>: :=  <single  space>|/|' 
<character>: :=  <letter> | <digit> | <special  character> | . 
<sign>::=  +|-|<empty> 

<identifier>:  :=  <letter>  |  <identifier><letter>  |  <identifierxdigit  > 
<character  code>::=  <character> | <character  code><character> 
<unsigned  integer>::=  <digit> | <unsigned  integer><digit> 
<integer>::=  <sign><unsigned  integer> 
<unsigned  number >: :=  <unsigned  integer > | .<unsigned  integer> 

|<unsigned  integer> .<unsigned  integer> 
<number>::=  <sign><unsigned  num"ber> 
<class  identifier>: :=  <identifier> 

<data  element  identifier>: :=  <identifier> |<class  identifier><identifier> 
<numeric  data  element  identifier>: :=  <data  element  identifier> 
<data  element  value>::=  <integer> I <number> I <character  code> I "<character  code>" 
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2.   Functional  Transformations 
<function  request>::=  <the  part>  FUNCTION  <function  identifier><is  part> 

<function  definition># 
<function  identifiers :=  <identifier> 
<function  definition>: :=  <coordinate  pair><coordinate  pair> 

| <coordinate  pair><function  definition> 
<coordinate  pair>::=   (<coordinate><limit  indicator >, <coordinate>  ) 
<coordinate>: :=   <number> 
<limit  indicator>: :=  +|-|<empty> 
<the  part>::=  THE|<empty> 
<is  part>::=  IS|<empty> 
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3.  Arithmetic  Expressions  Within  an  Occurrence 
<occurrence  arithmetic  expressions  :=  <occurrence  term> 

| < unary  operator><occurrence  term> 

| <occurrence  arithmetic  expression><add  operator><occurrence  term> 
<unary  operator>::=  +|- 
<add  operator>::=  +|- 
<occurrence  term>::=  <occurrence  factor > 

|<occurrence  term><multiply  operator ><occurrence  factor> 
<mult  iply  operator  > : : =  * | / 
<occurrence  factor>::=  <occurrence  primary> 

|<occurrence  factor><exponentiate  operator><occurrence  primary> 
<exponentiate  operator>::=  ** 
<occurrence  primary>::=  <unsigned  number> |<occurrence  variable> 

|<occurrence  function> | (<occurrence  arithmetic  expression>) 
<occurrence  variables  :  =  <numberic  data  element  identified 
occurrence  function  ::=  FUNCTION [<f unction  definition>] 

(<occurrence  arithmetic  expression>) 

|<function  identifier >  (<occurrence  arithmetic  expression>) 
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k .      Boolean  Expressions  Within  an  Occurrence 

occurrence  Boolean  expressions  :=   occurrence  Boolean  term> 

|<occurrence  Boolean  expression>OR<occurrence  Boolean  term> 
<occurrence  Boolean  term>::=  <occurrence  Boolean  factor> 

[<occurrence  Boolean  term>  AND  <occurrence  Boolean  factor> 
<occurrence  Boolean  factor>::=   <occurrence  relation> 

| (<occurrence  Boolean  expression>) 
<occurrence  relation>::=   <data  element  identif ier><equality  relational 
<data  element  value> 

|<data  element  identifier><one  of  relational>(<data  element  value  list>) 

|<occurrence  arithmetic  expression><relational> 
<occurrence  arithmetic  expression> 

|<occurrence  arithmetic  expression><relational><number> 

| <number><relational><occurrence  arithmetic  expression> 

| -Occurrence  arithmetic  expression><one  of  relational>(<integer  list>) 
<equality  relational>: :=   Is|lS  NOT|<is  part>  EQUAL  <to  part> 

|<is  part>  NOT  EQUAL  <to  part>|D0ES  NOT  EQUAL 

| EQUALS | EQ | NE | NQL | NEQ 
<is  part>: :=   is|<empty> 
<to  part>::=  TO|<empty> 

One  of  relational>:  :=   <is  part>  ONE  OF|<is  part>  NOT  ONE  OF 
<relational>: :=   <equality  relational > 

|  <is  part>NOT  <gtr  lss  part><than  part> 

|<is  part><gtr  lss  part><tha.n  part> 

|  <relational>  OR  <relational> 
<gtr  lss  part>::=  GREATER |GT |GTR | LESS |LT | LSS 
<than  part>::=  THAN|<empty> 
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<data  element  value  list>::=   <data  element  value> 

| <data  element  value>,<data  element  value  list> 
<integer  list>::=  <integer> | <integer>,<integer  list> 
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5.   Arithmetic  Expressions  Within  a_  Class 
<class  arithmetic  expression> : := 

<functional  designator><occurrence  arithmetic  expression> 
|<functional  designator><occurrence  arithmetic  expression> 
<occurrence  where  clause> 
<occurrence  where  clause>::=  WHERE  <occurrence  Boolean  expression> 
<functional  designator>: := 
SUM | TOTAL 
| AVE | AVERAGE 
| MIN | MINIMUM 
| MAX | MAXIMUM 
I  COUNT 
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6.   Arithmetic  Expressions  Over  a_  Parcel 
<parcel  arithmetic  expressions  :=   <parcel  term> 

| <unary  operator ><parc el  term> 

|<parcel  arithmetic  expression><add  operatorxparcel  term> 
<parcel  term>::=  <parcel  factor> 

|  <parcel  termxmultiply  operatorxparcel  factor> 
<parcel  factor>::=   <parcel  primary> 

|  <parcel  factor ><exponentiate  operatorxparcel  primary> 
<parcel  primary>::=  <unsigned  number > | <parcel  variable> 

|<parcel  function> | (<parcel  arithmetic  expression>) 
<parcel  variables  :  =   <class  arithmetic  expression> 
<parcel  function>::= 

FUNCTION  [<function  definition>] (<parcel  arithmetic  expression>) 

|<function  identifier>(<parcel  arithmetic  expression>) 
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7.   Boolean  Expressions  Over  a_  Parcel 

<parcel  Boolean  expressions :=  <parcel  Boolean  term> 

|<parcel  Boolean  expression>  OR  <parcel  Boolean  term> 
<parcel  Boolean  term>::=   <parcel  Boolean  factor> 

| <parcel  Boolean  term>  AND  <parcel  Boolean  factor> 
<parcel  Boolean  factor>::=   <occurrence  Boolean  expression> 

|(<parcel  Boolean  expression>) 
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8.   Regions 
<region  requests  :  = 

<the  part>  REGION<region  identifierxis  partxregion  expression:^ 
<region  identifier>: :=  <identifier> 
<region  expressions  :  = 

<region  specification> 

|  <region  specificationxset  operationxregion  specification> 

|  <region  specificationxparcel  where  clause> 
<parcel  where  clause>::=  WHERE<parcel  Boolean  expression> 
<set  operation>: := 

INTERSECT | OVERLAP 
| UNION | COMBINE 
| EXCLUDE | WITHOUT 
<region  specif ication>:  :=  <the  partxparcel  specification  list> 

|  <the  partxregion  partxregion  identifier> 
<region  part>::=  REGION | <empty> 
<region  identifiers  :=  <identifier> 
<parcel  specification  list>::=  <parcel  specification> 

| <parcel  specif ication>,<parcel  specification  list> 
<parcel  specif ication>:  :=  <qqsec  listxsectionsxtownships> 

| <sectionsxtownships> | <townships> 
<qqsec  list>::=   <qqsec  unit> 

| <qqsec  unit>,<qqsec  list> 

| <qqsec  unit>&<qqsec  list> 

|<qqsecs  unitxqqsec  list> 

| (<qqsec  list>) 
<qqsec  unit>::=  SEQ| SWQ | NEQ | NWQ | NH | SH|EH|WH 
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<sections>: :=   <section  numbers> 

SECTION  <section  numbers>|SEC  <section  numbers> 

<section  numbers>::=  <section  interval^ 

|<section  interval> ,<section  numbers> 
<section  interval>::=  <section  number> 

|<section  number>-<section  number> 
<townships>: := 

T<township  numbers xtownship  direction>R<range  number s><range  direction> 
<township  numbers>::=   <township  number> 

|<township  number >-<township  number> 
<township  directions  :=  N|s|u|L 
<range  numbers>::=   <range  number > 

| <range  number>-<range  number> 
<range  direction>: :=  e|w 
<section  number>::=  l|2|...|36 
<township  number > : : =  0 | 1 | 2 | . . . | 63 
<range  number>::=  0|l|2J...|63 
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9.   Primary  Requests 

<request>::=  <primary  request># | < secondary  request># 

<primary  request>::=  <on  phrasexfor  clause><keyword><object> 

| <on  phrase ><keyvord><f or  clause><object> 

| <on  phrase><keyvord><object><for  clause> 

|  <for  clause><on  phrase><keyword><o"bject> 

| <for  clause><keyvord><on  phrase><object> 

| <keyword><for  clause><on  phrase><object> 

|  <keyword><on  phrasexfor  clause><object> 

| <keyword><on  phrase><object><for  clause> 
<on  phrase>::=  ON  PRINTER 

|  ON  FILE  "<file  identifier>" | <empty> 
<file  identifier>: :=  <identifier> 
<keywor d> : : =  TABULATE | CALCULATE | MAP | OUTPUT 
<object>::=  <tabulate  object> 

| <calculate  object> | <map  object> | <output  object> 
<for  clause>::=  FOR<region  expression> | <empty> 
<tabulate  object>: :=  <tabulate  element > 

|<tabulate  element>  ,<tabulate  object>  |  <empty> 
<tabulate  element >: :=  <data  list> 

|  <data  listxparcel  where  clause> 
<data  list>::=   <data  specification> 

|  <data  specif ication> ,<data  list> 
<data  specif ication>: :=   <class  identified 

|  <data  element  identifier > 

I  occurrence  arithmetic  expression  > 
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<calculate  object>::=   <calc  partxwhere  phrasexby  phrase> 

|<calc  part><"by  phrase><where  phrase> 

| <where  phrase><calc  partxby  phrase> 

|  <where  phrasexby  phrase><calc  part> 

|  <by  phrase><calc  partxwhere  phrase> 

| <by  phrase><where  phrase><calc  part> 
<calc  part>: := 

<the  part><calc  list><of  part><occurrence  arithmetic  expression> 
<calc  list>::=   <calc  type>,<calc  list> 

|<calc  type>  AND  <calc  list> 
<calc  type>: := 
TOTAL | SUM 

|  AVERAGE  |  AVE 

| COUNT 

| PERCENTAGE | PERCENT | PCT 
<where  phrase>::=   <parcel  where  clause> | <empty> 
<by  phrase>::=  BY  <data  element  identifier> 

|BY  <occurrence  arithmetic  expression> | <empty> 
<map  object>::=   <parcel  arithmetic  expression> 

|AT  l:<unsigned  integerxparcel  arithmetic  expression> 

|<parcel  arithmetic  expression>  AT  l:<integer> 
<output  object>::=   <parcel  arithmetic  expression> 

|<parcel  arithmetic  expression^  <output  object> 


128 
10.   Secondary  Requests 
<secondary  requests>::=  <entity  remember  request> 

| <what  request > 

|<list  request> 

| <save  request > 

|<forget  request> 

| <make  request > 

|<leave  NARIS  request> 
<entity  remember  request>::=  <function  request> 

|<region  request> | <abbreviation  request> 
<abbreviation  request>::= 

<the  part>  ABBREVIATION  <abbreviation  identifier><is  part><expansion> 
abbreviation  identifier>: :=  <identifier> 
<expansion>: :=  <character> | <character><expansion> 
<what  request>::=  WEAT  <is  the  part><what  object> 
<is  the  part>::=  IS<the  part>  |ARE<the  part> 
<what  object >: :=   <entity  list> 

| VERSIONS | FUNDS | SPACE 

| <what  is  data  part> 
<entity  list>: :=   <entity  names  with  type> 

|<entity  names  with  type>,<entity  list> 

|<entity  names  with  type>  AND  <entity  list> 
<entity  names  with  type>::=  <entity  typexentity  name  list> 
<entity  type>::= 

REGION | REGIONS 

| ABBREVIATION | ABBREVIATIONS 

I  FUNCTION  I  FUNCTIONS 
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<entity  name  list>::=   <entity  narae> 

|<entit;y  name>,<entity  name  list> 

|<entity  name>  AND  <entity  name  list> 
<entity  name>::=   <region  identifier> 

abbreviation  identifier>  |  <function  identifier> 
<what  is  data  part>::=   <class  data  part> 
<class  data  part>::=   <data  part>  CLASSES 

| <class  identifier^ 

|<class  identifier><data  part><within  class  list> 
<data  part> : : =  DATA | <empty> 
<vithin  class  list>::=   <data  element  part> 

| <data  element  part>,<within  class  list> 

|<data  element  part>  AND  <within  class  list> 
<data  element  part>::=  ELEMENTS 

|<data  element  identifier^ 

|<data  element  identif ier><data  partxvalue  part> 
<value  part>::=  VALUES | <value  list> 
<value  list>: :=   <data  element  value> 

| <data  element  value>,<value  list> 

| <data  element  value>  AND  <value  list> 
<list  request>::=  LIST<the  partxlist  object> 
<list  object>::=  <drctry  entity  type> 

| <data  class  element  type> 
<drctry  entity  type>::=   <drctry  typexentity  type> 

| <drctry  type> | <entity  type> 
<drctry  type> : : =  PRIVATE | SEMIPUBLIC | PUBLIC 
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<data  class  element  type>::=  <data  part>  CLASSES 

|<class  identifier><data  part>  ELEMENTS 

j  <data  element  identifierxdata  part>  VALUES 
<save  request>::=  SAVE<entity  list> 
<forget  request>::=  F0RGET<entity  list> 
<make  request>::=  <drctry  type><entity  list> 
<leave  NARIS  request>::=  SUSPEND | SUSPEND<unsigned  integer>|END 
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